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1 .0  INTRODUCTION 


1.1  OBJECTIVE 


To  stucjy  active  structure  control  technology  as  applied  to  Large 
Space  Structures  (LSS).  In  particular,  to  develop  analysis  and  design 
techniques  for  utilization  of  practicing  control  system  engineers  In 
their  Investigations  of  future  Large  Space  Structures  (LSS).  These 
techniques  will  be  "use"  qualified  through  application  to  a  representative 
model  of  a  LSS.  This  objective  Is  in  support  of  the  overall  ACOSS  ob¬ 
jective  which  Is  to  develop  and  understand  a  generic,  unified  structural 
dynamics  and  control  technology  base  for  LSS  with  stringent  LIne-Of-SIght 
(LOS)  and  figure  performance  requirements  that  must  be  maintained  in  the 
presence  of  onboard  and  natural  disturbances.^ 


1.2  SCOPE _ 


^-=^Th1s  effort  will  be  restricted  to  theory  development,  experiment 
design,  demonstration  and  documentation  of  structural  dynamics  and  control 
technology  for  large  precision  spacecraft^  The  contractor  shall  plan  the 
overall  effort,  define  the  problem,  develop  and  evaluate  solutions.  The 
development  and  evaluation  of  solutions  will  Involve  analytical  studies, 
modeling  and  simulations.  The  contractor  shall  Identify  and  document  all 
the  analyses,  studies,  designs,  models,  simulations,  tradeoffs,  issues, 
software,  requirements,  specifications  and  data  produced  under  the  contract. 
Inportant  products  of  the  effort  will  be  the  demonstration  and  documentation 
of  new  control  technology. 


1.3  STATEMENT  OF  WORK 


The  following  four  tasks  comprise  the  original  Statement  of  Work.  Task  1 
has  been  deferred  for  implementation  in  AC0SS>17.  The  content  was  sub¬ 
sequently  modified  to  incorporate  Tasks  5  -  7  to  support  the  VCOSS  (Vibration 
Control  of  Space  Structures)  Program. 
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1.3.1  ORIGINAL  TASK/TECHNICAL  REQUIREMENTS 

Task  1  -  Investigation  of  Structural  Damping  Models:  Investigate  the  problem 
of  how  to  accurately  represent  structural  modal  damping  for  a  structure  that 
is  subjected  to  disturbance  and  control  forces.  Structures  with  values  of 
critical  damping  of  less  than  1%  can  be  adequately  described  by  assuming 
equivalent  viscous  modal  damping  coupled  with  real  eigenvectors.  However, 
for  larger  amounts  of  damping  these  assumptions  do  not  appear  to  be  true. 
Investigate  this  problem  in  depth  and  develop  a  complete  understanding  of 
the  proper  form  and  implementation  of  structural  damping  for  use  in  the 
analysis  and  development  of  active  control  technology  for  LSS. 

Develop  a  mathematical  model  containing  discrete  viscous  damping  at 
related  modes  throughout  the  structure. 

Compare  the  response  of  the  discrete  viscous  damping  model  to  that  of  an 
identical  model  containing  modal  rather  than  discrete  equivalent  viscous 
dampi ng . 

Develop  a  mathematical  model  with  an  accurate  form  of  structural  damping 
(proportional  to  displacement  and  at  a  90  degrees  phase  angle)  followed  by 
comparison  of  the  response  with  the  models  developed  above. 

Task  2  -  Development  of  Digital  Control  Analysis  Techniques;  Investigate 
sampling  phenomena  arising  from  a  digital  version  of  an  onboard  controller, 
using  a  truncated  version  of  the  model  developed  by  the  contractor  under 
Contract  F30602-81-C-0179.  Give  special  anphasis  to  the  effects  of  "folding" 
of  the  higher  neglected  modal  bending  frequencies  into  the  control  system 
bandwidth.  Establish  design  criteria  so  that  onboard  computer  sampling  rates 
can  be  selected  and  the  appropriate  tradeoffs  can  be  made. 

Task  3  -  Flexible  Appendage  Model:  Perform  and  document  modifications  of  the 
transfer  matrix  technique  for  simple  structures  with  rigid  appendages  that 
are  applicable  to  LSS.  Include  extensions  of  the  technique  that  contain 
flexible  appendages.  The  resulting  equations  shall  be  documented  with  their 
attendant  computer  program  implementations. 

Task  4  -  Stability  of  Large  Scale  Systems:  Use  the  decomposition  principle 
that  was  applied  under  Contract  F30602-81-C-01 79  to  stabilize  a  preliminary 
LSS  model  and  apply  it  toihe  higher  fidelity  models  that  are  being  developed 
under  Tasks  2  and  3. 
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Oral  Presentations:  Oral  presentations  shall  be  given  at  such  times  and 
places  as  designated  by  the  Contracting  Officer.  Sufficient  notices  will 
be  provided  prior  to  any  briefing. 
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1.3.2  "VCOSS"  TASKS/TECHNICAL  REQUIREMENTS 

Task  5  -  Develop  requirements  for  LSS  models;  Determine  system  characteristics. 
Assimilate  the  system  level  requirements  into  a  format  and  set  useable  by  a 
Structural  Analyst  to  provide  the  basis  for  an  accurate  and  reliable  structural 
model . 

Develop  derived  control  system  requirements.  Structure  the  system  level 
requirements  in  the  form  of  derived  requirements  as  viewed  by  a  Controls 
Analyst.  This  shall  include,  but  not  be  limited  to,  the  specification ;of 
potential  sensor  and  actuator  locations  and  the  limitations  oh  the  structural 
natural  frequency,  if  any. 

Determine  the  form  of  the  structural  models.  Review  the  requirements 
developed  above,  both  system  and  derived  control,  and  determine  the 
appropriate  technique  to  be  applied  to  the  models.  A  Lagrangian  approach 
shall  be  investigated.  It  may  not  be  used  if,  when  compared  with  other 
techniques,  it  proves  to  be  inferior  in  terms  of  taking  full  advantage 
of  modern  high  speed  computer  technology  and  maximizing  the  degree  of 
understanding  and  "feel"  of  the  expected  results. 

Task  6  -  Develop  "stiff"  structural  model;  Develop  mass  and  stiffness  matrices. 
Develop  mass  and  stiffness  matrices  for  the  candidate  structure  using  the 
techniques  and  procedure  in  Task  5.  Develop  a  model  with  sufficient  degrees 
of  freedom  and  fidelity  to  meet  contemporary  LSS  specifications. 

Develop  structural  characteristics.  Derive  the  structural  characteristics 
(eigenvalues  and  eigenvectors)  for  all  significant  modes  of  the  required 
model.  Develop  the  structural  gains  at  all  possible  sensor,  actuator,  and 
disturbance  locations.  If  multiple  paths  are  required  (as  in  multiple 
reflections  of  mirror  surfaces),  then  the  modal  gains  associated  with 
these  shall  also  be  developed. 

Task  7  -  Develop  "flexible"  structural  model:  Develop  mass  and  stiffness 
matrices.  Develop  the  mass  and  stiffness  matrices  for  the  candidate  structure 
using  the  technique  and  procedure  determined  under  Task  5.  Develop  a  model 
with  sufficient  degrees  of  freedom  and  fidelity  to  meet  the  contemporary  LSS 
specifications. 

Determine  eigenvalues  and  vectors.  Derive  the  structural  characteristics 
(eigenvalues  and  eigenvectors)  for  all  significant  mods  of  the  required  model. 
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Develop  the  structural  gains  at  all  possible  sensor,  actuator,  and  disturbance 
locations.  If  multiple  paths  are  required,  as  In  multiple  reflections  off 
mirror  surfaces,  then  the  modal  gains  associated  with  these  shall  be  developed. 


1.4  SCHEDULE 

The  original  scheduled  duration  of  the  contract  was  the  twelve  month  period 
beginning  15  April  1981.  As  a  result  of  VCOSS  modification,  the  VCOSS 
Final  Report  (Incorporated  within  this  Final  Report)  was  delivered  (with  an 
oral  report)  prior  to  30  December  1981.  The  ACOSS-15  Final  Presentation  was 
delivered  to  LTC  Herzberg  (DARPA)  and  Mr.  Carman  (RADC)  on  28  April  1982 
In  Huntsville,  Alabama.  A  compilation  of  due  dates  and  deliverable  Items 
comprises  Table  I. 


TABLE  I.  DUE  DATES  AND  STATUS  OF  DELIVERABLE  ITEMS 


DATE 

ITEM 

1  June  1981 

* 

R&D  Status  #1 

* 

ACOSS  15  Schedule 

★ 

ACOSS  15  Milestones/WBS 

.  ■  ■ 

1  July  1981 

* 

R&D  Status  §2 

1  August  1981 

* 

R&D  Status  #3 

1  September  1981 

* 

R&D  Status  #4 

1  October  1981 

* 

R&D  Status  #5 

’.-'i 

1  November  1981 

* 

R&D  Status  #6 

.  1 

•| 

1  December  1981 

* 

R&D  Status  #7 

-.  *•  • 

•  * 

1  January  1982 

* 

R&D  Status  #8 

* .  • 

1  February  1982 

* 

R&D  Status  #9 

1  March  1982 

* 

R&D  Status  #10 

1  April  1982 

* 

R&D  Status  #11 

1  Nay  1982 

* 

R&D  Status  #12 

* 

Final  Oral  Presentation 

* 

Final  Technical  Report 

* 

Abstract  of  New  Technology 

.1 
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1.5  CONTRACT  IMPLEMENTATION 

Dr.  Seltzer  was  the  Principal  Investigator  for  the  ACOSS-15.  He  also  perfonned 
the  Task  2. 

Dr.  Worley  and  Dr.  Glaese  perfonned  the  main  effort  performing  Tasks  3,  5*  6, 
and  7. 

Dr.  York  performed  the  main  effort  concerning  Task  4. 

The  effort  associated  with  Task  1  was  deferred  for  performance  under  ACOSS-17. 
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VCOSS  FLEXIBLE  APPENDAGE  MODELING 


Tasks  3  ("Flexible  Appendage  Model"),  5  ("Develop  requirements  for  LSS 
models"),  6  ("Develop  'stiff  structural  model"),  and  7  ('Develop  'flexible' 
structural  model")  are  reported  upon  in  one  cohesive  appendix.  This 
appendix  is  a  report  entitled  "DARPA-WPWFB  VCOSS  Dynamic  Model"  and  was 
presented  to  DARPA  and  Wright  Patterson  AFB  personnel  (as  well  as  Lockheed 
and  TRW  ACOSS  contractor  personnel)  on  1  October  1981. 

The  report  presents  analysis  results  applied  to  a  model  developed  for  DARPA 
by  the  Charles  Stark  Draper  Laboratory.  The  analysis  technique  embodies  a 
combination  of  finite  element  modeling,  closed-form  solutions  of  continuous 
beams,  and  modal  synthesis.  The  advantages  of  this  approach  are  that 
detailed  modeling  can  be  applied  where  necessary  while  still  retaining  the 
ability  to  make  modifications  and  changes  with  a  minimum  of  effort  and 
elapsed  time. 
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3.0 


DEVELOPMENT  OF  DIGITAL  CONTROL  ANALYSIS  TECHNIQUE 


The  purpose  of  this  section  Is  to  describe  an  Integrated  analytical  approach 
to  the  development  of  digital  control  systems.  This  philosophy  Includes  the 
dynamic  effects  of  sampling  phenonena  early  In  the  design  development.  Further 
It  attempts  to  enhance  computer  simulation  efficiency  by  Incorporating 
analysis  tools  to  prescribe  numerical  values  for  system  parameters  and  to 
predict  simulation  results.  The  main  body  of  this  section  was  presented 
at  the  IEEE  EASCON  In  Washington,  D.C.  on  November  18,  1981. 

The  design  philosophy  espoused  In  this  section  takes  advantage  of  high-fidelity, 
high-order  digital  simulation  tools  coupled  with  low  order  analysis  techniques 
that  permit  mathematical  tractablllty.  After  describing  that  philosophy,  the 
major  portion  of  the  section  Is  devoted  to  a  description  of  three  analytical 
digital  design  tools  and  how  they  can  be  used  In  concert.  These  tools  are: 

(1)  the  Systematic  Analytical  Method  (SW)  which  Is  a  method  for  obtaining 
mathematical  relationships  between  desired  system  Inputs  and  system  outputs; 

(2)  the  Cross-Multiplication  Method  which  provides  a  means  of  obtaining  the 
time-domain  response  from  the  system  closed-loop  transfer  function;  and  (3) 
the  Digital  Parameter  Space  Method  which  permits  one  to  examine  and  specify 
the  dynamics  of  a  system.  The  use  of  these  tools  together  Is  demonstrated  i 
In  a  simple  example. 

Since  the  preponderance  of  control  system  theory  has  been  developed  from  the 
standpoint  of  linear  and  continuous  systems,  the  design  analysis  process  usually 
starts  with  the  assumption  that  the  system  Is  continuous  and  not  digitally 
sampled.  Thus,  the  Intrinsic  quality  of  digital  computation,  namely  sampling, 
tends  to  be  Ingnored  until  the  last  design  step  at  which  point  one Xtyptcal.Ty 
says,  “Well,  we'll  simply  sample  'fast  enough'  so  that  It  doesn't  make 
any  difference".  However,  as  the  control  system  becomes  more  complex  and 
additional  requirements  arise,  the  capability  to  sample  "fast  enough"  may  not 
be  available.  In  addition,  the  designer  must  answer  the  question,  "What  Is 
fast  enough?" 

In  this  section  It  will  be  demonstrated  that  the  sampling  phenomena  can  be 
Included  In  the  early  stages  of  the  design  process  to  yield  valuable 
Information  regarding  not  only  control  system  design  requirements  but 
vehicle  confinguratlon  requirements  as  well. 

The  design  philosophy  may  be  described  graphically  In  Fig.  3-1. 

Each  block  of  the  figure  will  be  described  below. 
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ESTABLISH  REQUIREMENTS  FOR 


3.1 .  REQUIREMENTS  AND  SYSTEM  CONCEPT 


First,  the  requirements  for  the  digital  control  system  must  be  clearly  set 
forth.  This  trivial  and  Innocuous-appearing  step  Is  often  overlooked  or 
poorly  performed,  contributing  to  great  quantities  of  subsequent  pain 
(to  the  designers)  and  sometimes  poor  performance.  After  the  requirements 
have  been  specified,  an  overall  control  system  concept  Is  developed.  At 
this  point  In  time,  two  paralled  approaches  may  be  Initiated:  The  high-fidelity 
simulation  approach  and  the  conceptual  development  approach. 
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3.2. 


HIGH-FIDELITY  SIMULATION  APPROACH 


This  usually  begins  In  the  form  of  a  finite  element  model  of  the  overall 
system.  The  outputs  of  this  model  are  usually  In  the  form  of  eigenvalues 
and  eigenvectors.  TThe  real  part  of  each  eigenvalue  represents  the  square 
of  the  respective  modal  frequency  (In  units  of  rad/s),  and  the  corresponding 
eigenvector  represents  the  modal  displacenent  for  each  of  the  Included  degrees 
of  freedom  (normalized  to  a  unit  generalized  mass)J  The  outputs  of  the  finite 
element  model*  I.e.  the  modal  descriptions,  are  then  used  In  the  development 
of  a  high-fidelity,  high-order  computer  simulation.  The  simulation  can  be 
all-digital  or  can  begin  In  that  form  and  evolve  to  a  simulation  that 
Includes  hardware-1n-the-loop^. 
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3.3 


CONCEPTUAL  DEVELOPMENT  APPROACH 


In  this  approach,  the  system  concept  is  simplified  into  an  analytically 
tractable  model.  This  model  is  then  analyzed  to  gain  an  understanding  of 
its  dynamics  in  detail.  It  also  has  the  purpose  of  providing  nimerical 
values  for  system  parameters  (e.g.,  control  gains,  sampling  periods,  etc.) 
to  use  in  the  high-order  computer  simulation.  Once  these  initial  goals  have 
been  achieved,  the  fidelity  (and  complexity)  of  the  model  is  gradually 
increased.  Ultimately,  the  dynamic  performance  of  this  model  is  compared 
with  that  of  the  hijgh-order  simulation  to  predict  the  latter's  character 
and  hence  assist  in  the  "de-bugging"  the  program.  If  dynamic  agreement 
does  not  accrue,  then  the  order  of  the  simplified  model  is  increased 
further  until  agreement  is  achieved. 

Concurrent  with  the  development  of  the  simplified  analytical  model,  a 
simplified  structural  model  is  developed.  The  purpose  of  this  structural 
model  is  to  provide  the  simplified  analytical  model  developer  with  modal 
characteristics  in  a  timely  manner.  Of  course,  if  it  is  felt  that  modal 
characteristic  will  be  iforthcoming  from  ttie  finite  element  model  in  a 
timely  (with  respect  to  the  control  system  designer's  simplified  analytical 
model  needs)  manner,  then  development  of  the  simplified  structural  model  is 
not  needed.  (In  the  author's  design  experience,  this  timeliness  has  never 
been  demonstrated.) 

Once  the  simplified  models  have  been  c(xnpleted  intergrated,  they  are  used  in 
the  development  of  a  time-response  computer  model.  The  purpose  of  this  model 
is  to  verify  the  overall  flow  and  compatibility  of  the  control  system 
concept.  Since  time-invariant  coefficients  are  assumed  in  the  differential 
and  difference  equations  used  when  the  simplified  analytical  design  tools  are 
applied,  the  simulation  provides  a  valuable  check  of  the  validity  of  these 
assumptions  during  simulated  operation. 


3.4  INTEGRATION  OF  THE  TVO  APPROACHES 

The  overall  control  system  concept  is  the  point  of  initiation  of  the  two 
approaches,  which  are  developed  in  paralled.  As  numerical  values  are 
developed  for  system  parameters  in  the  "Conceptual  Development  Approach," 
they  are  integrated  into  the  "High-Fidelity  Simulation  Approach"  develop¬ 
ment.  Further,  the  dynamic  response  and  performance  determined  by  using 
the  simplified  models  are  then  compared  with  the  responses  of  the  high- 
order  simulation  (and  its  various  subordinate  portions)  to  assist  in 
"de-bugging"  the  simulation. 


3.5  SIMPLIFIED  DESIGN  TOOLS 


The  primary  design  tools  used  in  the  "Conceptual  Design  Approach"  are: 

*  The  Systematic  Analysis  Method  (SAM).  This  is  a  method  for  obtaining 
mathematical  relationships  between  desired  system  inputs  and  system 
outputs.  These  relationships  may  then  be  used  to  investigate  and  prescribe 
desired  transient  and  steady  state  system  performance  characteristics  (by 
applying,  for  instance,  the  Digital  Parameter  Space  Method  and  the  Final 
Value  Theorem).  They  may  also  be  used  to  determine  the  system  response  in  the 
time-domain  (by  applying,  for  instance,  the  Cross-Multiplication  Method).  The 
SAM  technique  is  described  below. 

*  The  Cross-Multiplication  Method.  This  is  a  technique  for  obtaining  the 
response  of  a  digitally-controlled  system  from  its  closed-loop  transfer 
function.  Since  this  technique  has  already  been  described  in  a  similar 
IEEE  conference,  it  will  only  be  described  herein.* 

*  The  Digital  Parameter  Space  Method.  This  technique  may  be  used  to  examine 
and  specify  the  dynamics  of  the  digitally- control  led  system.  Briefly  it 
permits  one  to  map  the  locations  of  the  roots  of  the  system  characteristic 
equation  into  a  space  whose  coordinates  are  the  system's  free  parameters. 

It  is  similar  to  a  multi -dimensional  root  locus  plot.  This  method  (or  any 
other  method  the  control  system  designer  chooses)  is  used  to  determine  the 
system  parameters  and  digital  computer  sampling  period.  Because  the  technique 
is  described  in  the  literature,  it  will  only  be  summarized  in  this  paper. * 
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3.5.1  SYSTEMATIC  ANALYTICAL  METHOD  (SAM) 

SAM  is  implemented  by  performing  the  following  four  steps.  If  the  equations 
resulting  from  the;first  three  steps  are  placed  in  a  table  of  three  columns 
(one  for  each  step)»  they  are  easily  manipulated  to  perform  the  fourth  and 
final  step. 

1.  Step  No.  1  Obtain  "System  Equations" 

First  one  selects  as  unknowns  the  variables  at  the  inputs  to  each  of  the 
samplers  in  the  system.  The  original  equai.:ons  describing  the  system  are 
rewritten  in  the  Laplace  domain  and  tabulated.  If  the  system  is  already 
described  by  a  block  diagram,  the  "system  equations"  may  be  written  upon 
inspection. 

2.  Step  No.  2  Obtain  "Modified  System  Equations" 

If  any  of  the  "system  equations"  contain  terms  that  consist  of  the  product 
of  any  unsampled  system  variable  and  an  unsampled  transfer  function,  they 
must  be  modified  by  substitution  to  eliminate  the  unsampled  variables.  For 
orderliness,  these  modified  equations  are  also  tabulated  (in  the  same  rows 
as  their  respective  original  equations). 

3.  Step  No.  3  Obtain  "Pulsed  System  Equations" 

Pulse  transforms  are  now  taken  of  each  side  of  the  "modified  system  equations," 
yielding  "pulsed  system  equations."  If  the  asterick  (*)  is  used  to  denote 
the  ideal  impulse  sampling  phenomenon,  then  the  transformation  into  the  sampled 
domain  is  performed  acccording  to  the  relationships, 

(R)  *  =  R*  (III-l) 

(RG)  *  =  (III-2) 

(RG*)  *  =  R*G*  (III-3) 

4.  Step  No.  4  Obtain  Desired  Input/Output  Relationships 

All  system  variables,  both  pulsed  and  continuous,  may  be  found  in  either  the 
"system  equations"  (Step  No.  1)  or  the  "pulsed  system  equations"  (Step  No.  3). 
The  desired  output(s)  may  be  solved  for  by  solving  the  appropriate  "system" 
or  "pulsed  system"  equations,  substituting  as  necessary. 

The  pulsed  (starred)  forms  of  the  variables  and  transfer  functions  may  all  be 
transformed  into  the  z-domain  by  noting  that  the  z-transform  of  the  variable 
(or  transfer  function)  C*{s)  is  merely  C(z),  using  standard  digital 
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nomenclature.'*  (Recall,  z=e'’‘,  where  s  is  the  complex  Laplace  operator  and  T 
represents  the  digital  sampling  period.)  If  it  is  desired  to  find  an  output 
expression  in  modified  z-transform  form  (so  that  the  system  response  ultimately 
can  be  determined  at  any  instant  of  ttme),  that  may  also  readily  be  accommodated 
Suppose  a  pulsed  variable,  C*,  is  equal  to  the  product  of  an  unstarred  quantity 
(such  as  a  transfer  function),  ACs),  and  a  starred  quantity  (such  as  an  input 
vartable),  B*(s).  f.e, 

where  A(s)  or  B*(s)  may  be  equal  to  unity.  The  modified  z-transform  may  al¬ 
ways  be  obtained  from  Eq.  (III-4)  by  performing  the  transformation: 

C(z,m)  =  A(z,m)  B(z)  (111-5) 

where  A(z,m)  represents  the  modified  z-transform  of  the  quantity  ACs),  and 
B  (z)  represents  the  ordinary  z-transform  of  the  quantity  BCs).**  This  technique 
appears  to  be  an  attractive  alternative  to  obtaining  modified  z-transforms 
through  Signal  Flow  Graph  techniques  (which  may  of  course  be  done). 


3.5. 2.4  SYSTEM  RESPONSE  BY  CROSS-MULTIPLICATION  METHOD ^ 

Once  the  relationship  between  the  desired  output  and  the  desired  Input  has  been 
established  (e.g.  by  Step  No.  4  of  SAM)  In  the  pulsed  (starred)  form  and  trans¬ 
formed  Into  the  z  or  modified  z-domain,  that  relationship  may  be  expressed  as  a 
ratio  of  two  polynomlnals  In  z  (or  modified  Zt  If  desired).  For  example.  If  the 
desired  output  Is  denoted  as  C(z}  and  the  desired  Input  Is  denoted  as  R(z), 
their  relationship  In  the  z-domaIn  may  be  written  as 


C  (z) 
R  (z) 


(III-6) 


where  coefficients  aj  and  b|(  represent  the  system  parameters.  As  demonstrated 
In  Reference  2,  the  procedure  for  finding  the  response,  C,  in  the  time  domain 
merely  consist  of  three  steps.  First  one  crossmultiplles  the  numerators  on  each 
side  of  Eq.  (III-6)  by  the  denominator.  Then  each  term  of  the  resulting 
expression  Is  divided  by  the  expression,  Finally,  the  Real  Translation 

Theorem**  Is  applied  to  that  result,  yielding  values  of  the  desired  response  at 
each  sampling  Instant  nT,  I.e. 

f  M  N  7 

c(nT)  =  1_  7  Z  a.r  (n  -  N  +  j)T  -  Z  b.c  (n  -  N  +  k)  T f  (III-7) 
bn  IjK)  J  k.O  ‘  1 

where  j,k,M,N,  and  n  are  integers. 

If  It  Is  desired  to  find  the  intra-sampling  response,  it  may  be  accomplished 
by  applying  either  the  Submultiple  Method  or  the  Modified  z-Transform  Method. 
Both  techniques  are  merely  modest  extensions  of  the  above-described  Cross- 
Multiplication  Method  and  are  described  In  Ref.  2. 
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3.5.3  DIGITAL  PARAMETER  SPACE  METHOD* 


Briefly,  the  method  permits  one  to  map  the  locations  of  the  roots  of  the 
systems 's  characteristic  equation  onto  a  plane  (or  three-space,  if  desired) 
whose  coordinates  are  the  system's  free  parameters.  In  addition  to  the 
characteristic  root  locations,  such  things  as  contours  of  constant  relative 
damping  ratios  (c)  and  specified  exponential  time  constant  may  be  portrayed 
as  contours  on  the  parameter  space.  Thus  in  a  manner  somewhat  reminiscent 
of  the  classical  root  locus  method,  a  portrait  may  be  presented  of  all 
pertinent  aspects  of  the  system's  transient  response  as  functions  of 
several  parameters  (rather  than  the  single  parameter,  open-loop  gain). 
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3.5.4  EXAMPLE 


A  simplified  model  of  a  typical  satellite  system  is  shown  in  Figure  3-2 
The  associated  system  equations  of  motion  may  be  written  in  vector-matrix 
notation  by  applying  SAM.  The  unknowns  are  selected  to  be  the  inputs  to  the 
samplers,  i.e.  e(t).  Applying  Step  No.  1  of  SAM,  the  equations  of  motion 
for  these  two  unknowns  plus  the  system  output,  x(t),  may  be  written  as: 

X*6pFGf,Q  U*  (III-8) 

U  =  D  E*  (III-9) 

E  =  R-  Y  =  R-  CX  (III-IO) 

Capital  letters  are  used  in  Eqs.  (II 1-8)  -  (111-16)  to  denote  matrices  and  vectors, 
and  asterisks  (throughout  the  paper)  represent  impulse  (i.e.  ideal  or  perfect) 
sampling.  These  equations  may  be  tabulated  as  the  first  column  of  the  Table. 

Next,  Step  No.  2  of  SAM  is  applied  to  Eqs.  (III-8)  -  (III-IO)  to  obtain  the 


"modified  system  equations",  i.e. 

X  =  Gp  F  U*  (i.e.,  no  change)  (III-8) 

U  =  D  E*  (no  change)  (III-9) 

E=R-Y=R-(C6pF  U  *  (III-IO) 

Equations  (III-8),  (III-9),  and  (III-IO)  may  be  placed  in  column  2  of  the  Table. 
Next  Step  No.  3  is  applied  to  these  equations  to  obtain  the  corresponding 
"pulsed  system  equations"  (which  may  be  placed  in  column  3  of  the  Table): 

X*  =  (Gp  G  Gj,q  )*  U*  (III-12) 

U*  =  0*  E*  (III-13) 

E*  =  R*  -  Y*  =  R*  -  (C  Gp  F  Gj^q)*U*  (III-14) 


Finally,  Step  No.  4  may  be  applied  to  obtain  the  relationship  between  the 
system  pulsed  output,  X*,  and  the  referenced  pulsed  input,  R*,  by  substituting 
Eq.  (III-14)  into  Eq.  (III-13)  and  then  substituting  the  resulting  equation 
into  Eq.  (III-12),  ultimately  yell  ding 

X*  -  GpF  Gf,g*  D*  R*  (I  +  0*  C  GpF  6,^^*)’*  (III-15) 

where  I  represents  the  identity  matrix.  The  z-transform  of  each  side  of 
Eq.  (III-15)  may  be  taken  to  yield 

X(2)  =  GpF  Gf,o(2)  D(z)  R(z)  [I  +  0(z)  ^^6^(2)]*^  (III-16) 


-i 
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A  digital  proportional  Integral  differential  (PIO)  controller  may  be  postulated. 
If  the  Integration  operation  Is  approximated  by  the  trapezoidal  approximation. 


Oj(z) 


T(zn) 

2(z-l) 


(III-17) 


The  PIO  controller  0(z)  may  be  written  as 


D(z) 


U(z) 

E(z) 


Kp  Ojlz)  +  K^0/(z)  + 


(III-18) 


V^KpT(z^-l)  lAK^T^(z-l)^  ♦  IC^|(z-l)^ 
(z-1)^ 


If  one  defines  Ap,  A^,  as 

Ap  =*  l^Kp  T^/J.  A^  •‘mi  T^/J.  Ap  «‘Kj  T/J  (III-19) 

one  may  readily  obtain  the  closed-loop  transfer  function  described  by  Eq. 
(III-20)  and  Table  3-1: 

2 

I  dj  z^ 

X(z)  ^  J^O _  (III-20) 

R(z)  *  3 

I  gj  2^ 
j-o 


Table  3-1;  CLTF  Coefficients 


The  closed-loop  transfer  function  of  Eq.  (111-20*)  can  now  be  used  to 
determine  both  the  steady-state  and  the  transient  behavior  of  the  system. 
Further,  the  Cross-Multiplication  Method  may  be  applied  to  Eq.  (III-20) 
to  determine  the  time-domain  response  of  the  system. 

The  form  of  the  controller  selected  can  be  checked.  In  a  coarse  sense,  by 
determining  the  system  steady  state  response  to  an  expected  form  of  disturbance. 
This  can  be  done  by  applying  the  Final  Value  Theorem  to  Eq.  (III-20)  after  sub¬ 
stituting  the  z-domain  version  of  the  reference  Input  Into  that  equation  and 
solving  for  0.  If  a  typical  reference  Input  Is  a  unit  step  function,  for 
example,  then 

6^(2)  •  2/(2-l).  (III-21) 

Substituting  Eq.  (III-21)  into  Eq.  (III-20),  solving  for  §(z),  and  finally 
applying  the  Final  Value  Theorem**  to  the  result,  one  obtains  the  steady  state 
value  for  0(t): 

11m  6(t)  -  lira  {(l-z’b(M2))*(fi(2)/0  (z))}  -1  (III-22) 

t-H*  z-1 

Not  unexpectedly,  this  Is  the  desired  value  In  the  steady  state. 

One  may  now  solve  for  the  time  response  of  the  system  by  applying  the  Cross- 
Multiplication  Method  to  Eq.  (I 11-20),  yielding 

((Ap+A^)ec(k-1)  +  2A^0g(k-2) 

(Ap-A.)0^(k-3)  -  (Ap+A^_A,,-3)0{k-l) 

(2A^-2Aj43)S(k-2)  +  (Aj,-A^-Aj+l)e(k- 

where  the  symbol  0(k)  represents  the  value  of  §(kT)  at  the  kth  value  of  T, 

I.e.  at  the  time  instant  k=kT.  Note  that  the  difference  equation  (III-23)  can 
admit  any  Initial  conditions  and  that  the  form  of  the  reference  Input, 
need  not  be  specified  except  at  the  sampling  Instants  (I.e.  it  can  be  complicated 
In  form  or  even  irrational). 

Finally,  the  nunerical  values  of  the  control  system  gains  and  the  sampling 
period,  T,  may  be  specified  with  the  aid  of  the  Parameter  Space  Method. 

Because  It  has  been  amply  described  In  the  literature,  only  a  final  form 
of  a  parameter  plane  Is  shown.  An  example  of  the  use  of  this  technique  Is 
briefly  shown.  The  denominator  of  the  closed-loop  transfer  function,  Eq. 

(III-20),  becomes  the  system  characteristic  equation  when  set  equal  to  zero. 
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The  stablll^  boundaries  from  the  z-domaln  are  mapped  onto  a  selected  par¬ 
ameter  space  in  Aq,  Ap,  Aj;  the  plane  of  the  paper  defines  the  Aq  -  Ap  space 
(see  Figure  111-3).  The  z*+l  boundary  Is  defined  by  the  plane  Ap»0,  and 
the  z*-l  boundary  by  the  plane  Ao*2.  The  z®e^'^n^  boundary  Is  a  surface, 
the  contours  of  which  are  defined  by  curves  designated  by  their  associated 
nunberical  values  of  Aj  (as  Aj  varies  from  0  to  1.4). 

Contours  associated  with  various  values  of  damping  ratio  (()  are  then  mapped  onto 
the  same  space  (Figure  I I 1-4),  as  are  lines  associated  with  various  values  of 
real  root  locations  (6).  Since  the  characteristic  equation  Is  third  order,  one 
would  expect  the  space  to  be  spanned  by  Intersections  of  contours  Indicating 
locations  of  three  roots.  For  example,  at  the  suggested  design  point,  one  pair 
of  complex  conjugate  roots  has  a  damping  ratio  of  0.707  and  the  real  root  Is 
located  at  approximately  0.1.  From  this  frequency  domain  Information  the  system 
dynamics  can  be  deduced  In  terms  of  the  system  gains  (parameters)  and  sampling 
period  (T).  The  actual  time  domain  response  Is  then  found  by  using  the  Cross- 
Multiplication  Method  to  produce  and  appropriate  difference  Eq.  (111-23),  In¬ 
serting  the  numerical  values  associated  with  the  candidate  design  point 
selected  on  the  parameter  space.  Using  this  equation,  the  system  response  to 
various  forms  of  reference  Inputs  can  then  be  determined  explicitly  for  various 
choices  of  gains  of  sampling  period. 
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PARAMETER  PLANE  (AI=0.2) 


PARAMETER  PLANE  (Aj  *  0-2) 

FIGURE  III-4 

CANDIDATE  DESIGN  PT:  Ap  =  0-87,  Aj  =  0-2,  Aq  =  1.66 

«„T  =  1.93,  c  =  0.707,  fi  =  0.1 
X(T)  =  Ap  +  Aj  +  Aq  =  2.73 
X(2T)  =  0.55 
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STABILITY  OF  LARGE  SCALE  SYSTEMS 


4.1  INTRODUCTION 

The  problem  attacked  In  this  section  Is  to  develop  a  method  to  determine  stabi¬ 
lity  regions  of  feedback  control  gains  for  a  large  scale  OARPA  type  system,  such 
as  that  described  In  Section  2.0.  The  approach  Is  to  extend  the  parameter  space 
technique  for  discrete  control  systems  to  permit  Inclusion  of  an  arbitrary  number 
of  bending  modes.  The  large  scale  system  will  be  represented  In  modal  coor¬ 
dinates,  and  the  discrete  feedback  controller  will  employ  both  attitude  (6)  and 
attitude  rate  (S).  The  problem  will  be  to  find  the  respective  control  gains  Kp 
and  (see  Figure  4-1)  that  will  Insure  stability.  The  resulting  algorithm  to 
yield  the  required  form  of  the  system  characteristic  equation  must  be  suitable 
for  computer  Implementation  and  simulation. 

It  Is  assumed  that  the  digital  onboard  controller  will  be  of  the 
position/derivative  of  position  (PO)  form.  The  two  parameters  selected  for 
portrayal  on  the  parameter  plane  are  Kp  and  K^.  For  the  simplified  example,  only 
one  sensor  and  one  torquer  are  used,  and  colocation  Is  not  assumed. 


4.2  PARAMETER  SPACE  METHOD 

Based  primarily  on  the  work  of  D.  D.  Slljak,^  the  parameter  space  method  Is  an 
analytical  tool  developed  for  use  In  control  system  analysis  and  synthesis. 
Briefly,  the  method  as  developed  by  S.  M.  Seltzer,^  permits  one  to  map  the  loca¬ 
tion  of  the  roots  of  the  system's  characteristic  equation  onto  a  plane  (or  three- 
space  If  desired)  whose  coordinates  are  the  system's  free  parameters.  The 
characteristic  equation  must  be  of  the  form, 

M 

C.  E.  »  I  [A^  ko  +  +  F^l  (4-1) 

1  =  1 


The  plane  Is  readily  divided  Into  regions  Identified  with  system  stability  and 
Instability.  The  method  can  then  be  used  to  provide  a  graphical  display  of  sta¬ 
bility  regions  for  two  specified  parameters,  such  as  kQ  and  k^.  The  free  parame¬ 
ters  need  not  be  gains  but  just  as  well  could  be  other  system  parameters.  In 
addition  to  the  characteristic  root  locations,  other  factors  of  Interest 
such  as  constant  damping  ratios  and  specified  exponential  time  constant  may 


be  portrayed  as  contours  on  the  parameter  space.  Thus  In  a  manner  somewhat  rem1> 
niscent  of  the  classical  root  locus  method,  a  portrait  may  be  presented  of  all 
pertinent  aspects  of  the  system's  transient  response  as  functions  of  several 
parameters  (rather  than  the  single  parameter,  open-loop  gain),  and  as  a  function 
of  the  Independent  argument,  uT,  where  u  Is  the  system  damped  frequency  and  T  Is 
the  digital  sampling  period.  An  automated  generation  of  the  arrays  A,  B,  and  F 
for  an  arbitrary  number  (M)  of  bending  modes  Is  highly  desirable,  where 

A  »  {A^} , 

(4-2) 

F  *  {F^}. 

4.3  MODAL  TRANSFORMATION 

The  original  system  may  be  described  by  the  vector-matrix  equation. 


[m]  i  +  [k]  ^  =  F 


(4-3) 


where  [m]  Is  the  mass  matrix;  [k]  Is  ti'e  stiffness  matrix,  £  Is  the  force/ torque 
vector,  and  X_ Is  the  system  state  vector. 

Equation  (4-3)  may  be  transformed  Into  modal  coordinates  (n)  using 

^  =  q  T^,  (4-4) 

with  Q  the  modal  matrix  composed  of  the  normalized  eigenvectors. 

Substitution  of  Equation  (4-4)  Into  (4-3)  yields 

h  +  [u^]n  *  £=  Tg.  (4-5) 

The  torque,  T^.,  Is  of  the  form 


’’’c  *  '^p  ®  +  l^d  ®  • 
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(4-6) 


I  I  tl 


p 


h 


where  0  and  e  represent  angular  and  angular  rate  displacements.  It  Is  assumed 
that  the  rate,  0,  Is  measured  directly  and  that 


0  -  I 


♦j  '’j* 


(4-7) 


4.4  GENERAL  FORM  OF  CONTROLLER 

The  large  scale  system  represented  by  Equation  (4-3)  or  (4-5)  Is  to  be  controlled 
by  a  conventional  rate  (0)  plus  position  (o)  controller.  Knowledge  of  the 
damping  ratio  and  frequency  of  each  of  the  modes  Is  assumed.  The  general 
form  of  the  PO  controller  used  Is  shown  In  Figure  4-1. 

DIGITAL  COMPUTER 


FIGURE  4-1.  PROPORTIONAL  (0)  -  DERIVATIVE  (0)  CONTROL  LAWS 
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4.5  PLANT  DESCRIPTION 

Each  path  shows  the  transfer  function,  nj/T^. 

of  the  plant,  e/T^.  Is  given  by: 


and  the  overall 


transfer  function 


N 

e(s)  *  I  ♦j 

Tc(s)  Tc(s) 


(4-8] 


This  plant  is  portrayed  graphically  in  Figure  4-2. 


4.6  DERIVATION  OF  THE  CHARACTERISTIC  EQUATION 

Using  standard  feedback  control  simplifying  techniques,  one  can  show  that  the 
transfer  function  corresponding  to  the  control  system  shown  In  Figure  4-1,  Is 
given  by 


9* 


K  R  * 


(4-9) 


1  +  [Kd  +  KpO;(z)]Gp^ 


where 


Gp{s)  =  I  ♦^Y^P^(s), 


(4-10) 


1  =  1 


Gp*(z)  =  Z[Gj,o(s)  Gp(s)] 


(4-11) 


P,i(s)  *  2  -  2 

^  S  +  ZC^U^S  + 


(4-12) 


D;{z)  = 


T(zn)  , 
2(z-l) 


(4-13) 


1  -  e 


-Ts 


(4-14) 


with  Y^,  4^  denoting  the  slope  of  the  1th  mode  at  torquer,  sensor  (respectively), 
and  denoting  the  damping  ratio,  natural  frequency  of  the  ith  mode.  The 

sample  time  Is  T. 
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Notation:  The  following,  frequently  appearing  quantities  In  the  derivation  are 
represented  by  the  g1 ven  symbol : 

0.^ 

2  COS  3^T  (4-15) 

s^  2  sin  3^T 


Before  deriving  the  characteristic  equation,  ttiree  lemmas  will  be  stated  which 
simplify  the  derivation.** 


Lemma  1. 


Z[Gho(s)  G(s)]  *  Z[(l  -  e-Ts)  G(s)/s] 

*  (1  -  z“M‘Z[  G(s)/s] 


(4-16) 


Lemma  2. 


Z[ 


_ b 

(s+a)^ 


z  e“^^s1n  b  T _ 

z  -  2ze"*^cos  b  T  + 


(4-17) 
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Lemma  3.  If  u  4:  0,  then 


s^+  2cws  +  ^  <0(1  • 

_ z  e“^“^»  sIn[(d(l-c^)^T] _ 

z  -  2ze"^“^cos[<i)(l-c^)^T]  +  e"^^“ 


(4-18) 


Proof: 

(First  complete  the  square  of  the  denominator,  then  combine  the  linearity  of  the 
Z  transform  with  Lemma  2  to  obtain  the  stated  result). 

Restating  Lemma  3  using  the  notation  In  Equations  (4-15)  yields 


Z[ 


s^+2cws+«^ 


.  es2 _ 

e  -  2ecz  +  e^ 


Finally,  If  <0  =  0  ,  then 
Lemma  4. 


Z[ 


Tz 

(z-1)^ 


(4-19) 


(4-20) 


Note:  Lemma  4  can  be  shown  to  be  the  limiting  case  of  Lemma  3  as  u-*-0  using 
L' Hospital's  Rule. 

With  the  aid  of  the  above  limits,  the  characteristic  equation  can  now  be  derived. 
Equation  (4-13)  will  replace  Dj-(z),  but  the  expression  for  Gp*(z)  will  have  to  be 
found . 
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Gp*(2)  =  Z[Gho(s)  Gp(s)] 

*  (l-z"M  Z[  Gp(s)/s  ]  ,  by  Lemma  1  (4-21) 

N 

*  il'Z'h  I  4^  Z[l/(s^+  25^w^s  +  «^^)] 

1=1 


[  h  Yi 


N 

Jz  ^  I  *1 

(z-l)^  1*2  ^1  z^  -  2e^CfZ  + 


using  Lemma  4  and  3. 


Next,  substitute  0(z),  Gp*(z)  Into  the  denominator  of  the  transfer  function 
0*/r*,  I.e.,  Equation  (4-9),  to  obtain  the  characteristic  equation. 


1  +  [Kd  +  Kp 


[♦l  Yi 


Tz 

(z-1)^ 


N 

+  I  *1  ^1  *  ei^iz 

1*2  ®1  z^  -  2e|C^z  +  e^^ 


(4-22) 


Cancelling  and  multiplying  by  z-1  yields 

1  +  [Kj(z-l)  .  Kp  J  (IH)]  •[  IHLI 
‘  (z-1) 


+ 


N 


I 

1*2 


llll 

^1 


ei^i _ 

z^  -  2e^c^  z  +  e^^ 


1*0. 
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Finally,  Me  obtain  the  C.E.  in  polynomial  form  by  combining  the  finite  sum 

inside  brackets  into  a  single  fraction  and  multiplying  through  by  the  denomina¬ 
tor.  Let 


H,  -  llim.  1.2.3.  ...M. 
^i 


then 


HI  H,. 

0*  (z-1)^  i=2 


X  -  2e^c^2  +  e^^ 


Combining  yields 

N 

m  *  [♦i^iT  •  n  (z^  -  2e^c^z  +  e^^) 

D*  1=2 

N  N 

+  I  H|^(z-1)^  •  n  (z^  -  2e^c^z  +  e^^lJ/D* 
k«2  i»2 

i+k 


(4-24) 


(4-25) 


(4-26) 


where 


N 

D*  «  (z-1)^  n  (z^  -  2e^c^z  +  e^^) 
1=2 


(4-27) 


Substituting  and  multiplying  through  by  D*  gives  the  C.E.  in  polynomial  form 
(albeit  not  the  most  convenient  form): 

D*  +  [Kjjiz-D  l^p  *  ^  (z+1)]n*  »  0.  (4, 
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or 


N 

(z-1)^  n  (z^  -  2e^c<z  +  ♦  [Kd(z-l)  +  Kp  j  (z+1)]  . 

i»2 

N 

[♦l^iT  n  (z^  -  2e^c^z  +  (4-29) 

1*2 


N  N 

+  I  (z-l)^  n  (z^  -  2e^c^z  +  e^^)]  »  0. 

k*2  1»2 

1+k 


4,7  ALGORITHMIC  OEVELOPMEMT  OF  PARAMETER  PLANE  FORM  OF  C.E. 

For  the  parameter  plane  analysis,  let  Icq,  be  the  two  parameter  of  Interest. 
Then,  the  characteristic  equation  must  be  In  tiie  form 

M 

C.E.  *  I  (A^ko  +  B^ki  +  F^)z^  .  (4-1) 

1*0 


For  this  application, 

ko  *  ,  ki  *  Kp  .  T/2.  (4-30) 

The  objective  then  will  be  to  develop  an  algorithmic  approach  that  will  trans¬ 
form  the  characteristic  equation  given  In  Equation  (4-15)  Into  the  form  of 
Equation  (4-1),  that  Is,  generate  the  arrays  A,  B,  and  F. 

The  following  algorithm  for  multiplying  polynomials  will  be  Instrumental  In 
achieving  the  above  goal.  As  can  be  seen  from  Equation  (4-29),  the  polynomial 
product 

N 

p2n(z)  “  n  (z^  +  a^z  +  b^)  (4-31) 

1*1 
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occurs  quite  frequently.  Suppose  that  the  resulting  polynomial  Is  of  degree  m. 
Let  p^  *  denote  the  coefficient  of  the  term  In  this  polynomial.  For  example. 
If  then 

P2(z)  “  P2,2Z^  +  P2.1*  +  P2.0 

where 

P2.2  ■  1*  P2.I  "  »!•  1^.0  '  '>1  • 

ALGORITHM  P;  Let  N  ^  2  be  specified.  Then  the  coefficients  of  the  polynomial, 

N 

^2N  ®  (z^  +  a^  z  +  1)^) 

1=1 


can  be  determined  recursively  by 
P2m,0  “  P2m-2,0  ‘^m 
P2m,l  *  P2m-2,1  ‘’m  *  P2m-2,0  «m 

P2m,1  *  P2m-2,1  •^m  P2m-2,1-1  *m  *  P2m-2,1-2 

for  1  =  2,3,  ...,2m-2. 

P2m,2m-1  *  P2m-2,2m-2  ^  P2m-2,2m-3 

P2m,2m  “  P2m-2,2m-2 
where  m  =  2,  3,  .. .M 

Notes:  1.  All  Intermediate  products  are  available  as  well  as  the  final  polyno¬ 
mial  -  the  order  of  calculation  Is 

P2(z)  given,  the  P4(z),  Pgiz) . 
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2.  Although  It  is  possible  to  derive  closed  form  expressions  for  the 
coefficients  of  P2ii|(z),  the  above  recursive  approach  is  more 
accurate  computationally. 


Proof: 


The  proof  is  a  standard  application  of  the  method  of  mathematical 
induction. 

The  desired  form  of  the  characteristic  equation.  Equation  (4-1),  can  now  be 
obtained  from  Equation  (4-29).  The  goal  will  be  to  generate  polynomials  D(z) 
and  H(z)  so  that  the  characteristic  polynomial  can  be  expressed  as 


0(z)  +  [Kjj(z-l)  +  K  -  (z+1)]  H(z). 


(4-35) 


The  general  procedure  is  outlined  in  Figure  4-3.  Note  that  D(z)  and  H(z)  have 
degree  2N  and  2N-2,  respectively. 

Each  step  in  the  procedure  is  discussed  with  computational  efficiency  in  mind. 
Step  1.  The  polynomial 


0(z)  5  D2n(z) 


(4-36) 


can  be  obtained  using  ALGORITHM  P  with  initial  values  of 


d2  =  1,  dj^  *  -2,  dQ  =  1, 


(4-37) 


Notati on i  Let  c^  denote  the  coefficient  of  z^  in  polynomial  C(z),  i.e.  for  D(z) 
of  degree  2N,  obtained  in  Step  1, 


D(z)  *  d2fj  z^**  +  dgfj.j  z2N”l  +  ...+d2Z^+djZ  +  dQ 


(4-38) 


Step  2.  Obtain  the  coefficients  of  G(z)  from  those  of  the  previously  calculated 
0(z).  Since 


0(z)  »  (z^  -  2z  +  1)  •  G(z), 


(4-39) 
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FIGURE  4-3.  PROCEDURE  FOR  OBTAINING  ARRAYS  A,  B,  AND  F. 


multiply  out  the  right-hand  side,  equate  coefficients,  and  solve  for  the  coef¬ 
ficients  Qj,  1*0,  1,  2,  2N-2  to  obtain 

92N-2  “  ‘*2N 

92N-3  *  *  ^32H-2 

92N-2-1  “  ‘^N-l  *  “  92N-1 


for  1  =  2,  3,  ...,  2N-2. 


Step  3.  Express  the  coefficients  of  In  terms  of  those  of  Diz),  K  given. 

To  do  this,  let 


N 

D(z)  *  (z-1)^  n  (z^  +  a^  z  +  b^) 
1*2 


N 

(z^  +  aj^  z  +  bjj)  [(z-1)^  n  (z^  +  a^  z  +  b^)] 

1*2 

1+k 


(4-41) 


(z^  +  ajj  z  +  bj^)  •  E^*^^(z) 


vrhere  (referring  to  Equation  (4-15)), 


a^  *  -  2e^Cf,  b^  *  e^^ 


(4-42) 


As  done  above  In  Step  2,  multiply  out  Equation  (4-42),  equate  coefficients,  and 
solve  for  coefficient  e^^*^^(z),  1*0,  1,  2,  ...,  2N-2  to  obtain 


“  <*2N 


®^'^^2N-3  *  ‘*2N-1  “  *k  ®^*^^2N-2 

®^*^^2N-2-1  “  <*2N-i  ■  h  ®^'^^2M-f  "  \  ®^*^^2M-l-i 
for  i-2,  3 . 2N-2. 

STEP  4.  Adding  correspondl ng  terms  yields  for  H(z), 

N 

^1  “  *1^1^91  ■*■  ^  (4-44) 

k»2 

for  1-2M-2,  2M-3 . 2.  1,  0. 


Step  5.  With  H{2)  and  0{z)  obtained,  the  characteristic  polynomial  becomes 

0(2)  +  [Kj  (2-1)  +  KpX  (z+l)]  ♦  H(2)  .  (4-45) 

2 

With  kQ,  kj  being  the  two  parameters  of  Interest  as  given  In  Equation  (4-30), 
the  arrays  A,  B,  and  F  can  be  obtained.  For  example,  the  coefficient  of  Is 
given  by 


d2  +  Kj  •  +  Kp  •  X  ^  ’^p  1  *^2  • 

2  2 

Collecting  the  appropriate  terms  as  shown  In  Equation  (4-1)  yields 

^2  *  ^1  “  ^2»  ®2  *  I'l  h2»  ^2  *  ^  (4-47) 

The  complete  arrays  are  given  In  the  following  Table  4.1. 
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Al 

Fi 

z2N 

0 

0 

<*2N 

z2N-l 

h2N-2 

h2N-2 

<‘2N-1 

z2N-2 

h2N-3  -h2N-2 

h2N-3  +h2M-2 

<12(1-2 

z2N-3 

• 

• 

h2N-4  -h2n-3 

• 

• 

h2N-4  +h2M-3 

• 

» 

<l2n-3 

• 

• 

• 

• 

z2 

• 

« 

hi  -  h2 

• 

• 

hi  +  h2 

• 

» 

d2 

zl 

ho  -  hi 

ho  +  hi 

<ll 

zO 

-ho 

+ho 

<io 

TABLE  4.1.  TWO  PARAMETER  FORM  OF  CHARACTERISTIC  EQUATION 

(ARRAYS  A,  B,  AND  F) 
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9 


* 


C 


1 

i ; 

I , 

I 

« 

I 


4.8  VERIFICATION 

In  previous  contractual  work  on  large  space  structures,  the  parameter  space 
approach  was  employed  with  one  bending  mode.  Subsequently,  a  second  bending 
mode  was  added.  The  results  obtained  with  the  new  algorithm  that  Is  designed  to 
handle  an  arbitrary  number  of  bending  modes  are  In  agreement  with  those  of  the 
previous  work. 

4.9  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  STUDY 

The  development  of  this  algorithm  removes  a  serious  drawback  In  the  application 
of  the  parameter  plane  approach  to  the  problem  of  discrete  control  of  large 
space  structures.  If  the  number  of  desired  bending  modes  exceeds  two,  then  the 
task  of  obtaining  the  necessary  form  of  the  characteristic  equation  becomes 
extremely  complex.  Once  the  modal  representation  Is  available,  the  automation 
of  the  generation  of  the  needed  arrays  A,  B,  F  eliminates  this  drawback. 

Previous  work  by  Asner  and  Seltzer  had  extended  the  parameter  plane  approach  to 
arbitrarily,  high  order  continuous  control  systems.  It  Is  anticipated  that  this 
algorithm  will  now  be  employed  in  conjunction  with  the  modal  representation  of 
the  DARPA  flexible  large  scale  space  structure  developed  by  Control  Dynamics. 

Several  areas  of  possible  extension  of  this  work  follow: 

1)  Include  sensor,  torquer  dynamics 

2)  Allow  for  more  than  one  sensor,  torquer  pair. 

3)  Include  structural  damping  at  the  start,  rather  than  Introducing  It 
later  Into  the  modal  representation. 

4)  Introduce  a  second  Integrator  Into  the  feedback  loop  with  a 
corresponding  control  gain.  Preliminary  work  Indicates  that  the 
same  line  of  development  can  be  used  and  that  the  needed  arrays  A, 
j,  C,  F  can  be  generated  by  the  computer. 


.  3 


1 


< 
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5.0 


CONCLUSIONS 


With  the  exception  of  Task  1  (Investigation  of  Structural  Damping  Models), 
all  remaining  six  tasks  were  completed  on  time  and  within  schedule.  The 
development  and  prescription  of  an  orderly  approach  to  the  implementation 
of  the  developed  digital  control  analysis  techniques  was  performed  under 
Task  2  and  reported  upon  in  Section  3.0.  Closely  associated  with  the  work 
under  Task  2  was  the  extension  (to  the  digital  domain)  of  the  decomposition 
principal  (Parameter  Space).  This  was  reported  upon  in  Section  4.0.  All 
modeling  tasks  (Tasks  3,  5-76)  are  reported  upon  in  Section  2.0.  The  unique 
analysis  technique,  embodying  a  combination  of  finite  element  modeling, 
closed-form  solutions  of  continuous  beams,  and  modal  synthesis,  is  reported 
upon  in  Section  2.0.  This  work  (Tasks  3,  5-7)  presents  a  means  of  better 
understanding  the  actual  dynamics  of  the  resulting  mathematical  model. 


ABSTRACT 


This  report  presents  analysis  results  of  future  Large  Soace  Structures  (LSS) 
that  can  be  used  by  VCOSS  Contractors  in  their  control  system  analyses.  The 
data  is  presented  in  the  form  of  modal  data  for  two  different  "configurations" 
of  a  candidate  LSS.  The  analysis  technique  utilized  is  a  combination  of  finite 
element,  uniform  beam  theory  and  modal  synthesis.  The  modal  was  developed  to 
be  not  only  representative  of  LSS  characteristics  but  also  to  be  readily  amenable 
to  modification  and  adjustment  as  required.  The  model  used  in  this  analysis  was 
derived  from  a  model  presented  in  DARPA  Report  #R-1404  "Active  Control  of  Space 
Structures  Interim  Report"  by  the  Charles  Stark  Draper  Laboratory,  Inc. 
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1.0  INTRODUCTION 

The  purpose  of  this  report  is  to  present  the  results  of  an  analysis  activity 
undertaken  to  develop  detailed  structural  models  of  two  candidate  LSS  models 
labeled  Stiffness  and  Strength.  The  analysis  procedure  is  presented  in  an 
overview  format.  A  later  technical  report  is  planned  to  cover  in  detail  the 
technical  drivation.  Results  are  presented  in  the  form  of  modal  deflections 
for  the  optical  structure  and  the  equipment  section  independently,  in  addition 
to  coupled  mode  shapes  for  the  complete  structure.  These  uncoupled  data  are 
presented  in  this  manner  to  facilitate  analysis  by  the  VCOSS  contractors. 


2.0 


LSS  QUALITIES/REQUIREMENTS 


What  is  a  Large  Space  Structure?  A  question  that  appears  to  have  an  obvious 
answer  --  A  Targe  structure  in  space.  However,  a  more  careful  examination  re¬ 
veals  that  the  answer  is  too  simple  and  does  not  get  to  the  heart  of  the  matter. 
The  space  program  has  spawned  a  number  of  "Large  Space  Structures"  in  the  past 
if  size  is  the  only  criteria.  The  Saturn  class  of  launch  vehicles  weighed  over 
5  million  pounds,  and  the  Skylab  weighed  almost  two  hundred  thousand  (200,000) 
pounds.  Both  programs  were  very  successful  in  the  1960 's  and  early  1970's.  Thus 
mere  physical  size  is  not  enough  to  completely  describe  LSS  even  though  it  must 
be  an  ingredient.  What  then  are  the  appropriate  characteristics?  It  is  out 
opinion  that  an  LSS  is  characterized  as  follows: 

SUMMATION  OF: 

LARGE  DIMENSION 
LOW  STIFFNESS 

PRECISION  POINTING  REQUIREMENTS 
HIGH  MODAL  DENSITY 
LOW  FUNDAMENTAL  FREQUENCY 
distributed  actuators  and  SENSORS 

However,  to  this  list  a  very  important  ingredient  must  be  added.  It  is  important 
to  realize  that  we  are  trying  to  extend  the  caoability  of  control  systems  to  as 
high  a  degree  as  possible.  Therefore  the  characteristic  that  makes  an  LSS  unique 
is  the  fact  that  it  includes,  in  addition  to  the  above, 

STATE-OF-THE-ART  CONTROL  SYSTEM  TECHNOLOGY  ! 

It  is  through  this  added  ingredient  that  the  full  dimension  of  LSS  is  obtained. 

It  is,  therefore,  incumbent  on  the  structural  analysis  to  develop  a  model  that 
not  only  has  large  dimension,  etc.,  but  also  has  the  ingredients  that  cause  a 
control  system  designer  difficult.  The  purpose  of  this  report  is  to  present  the 
results  of  an  analysis  activity  undertaken  to  accomplish  that  task. 
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We  have  undertaken  an  analysis  of  this  particular  structure  with  the  viewpoint 
of  developing  a  model  of  a  large  space  structure  that  would  be  of  sufficient 
detail  and  complexity  to  challenge  the  capability  of  control  systems  designers 
while  retaining  management  flexibility.  The  techniques  utilized  embody  a  com¬ 
bination  of  finite  element  modeling,  closed  form  solution  of  continuous  beams 
and  modal  synthesis.  Since  the  light  path  is  of  crucial  concern  in  this  analysis 
we  have  developed  an  accurate  model  of  the  motion  of  the  central  ray  as  a  function 
of  the  motion  of  each  of  the  elements  in  the  optical  path.  In  addition,  we  have 
constructed  our  model  so  that  the  isolation  section  stiffness  can  be  readily  mod¬ 
ified  as  requirements  dictate. 

Our  analysis  procedure  will  be  to  employ  modal  synthesis.  This  technique  has 
been  employed  successfully  in  the  past  on  a  number  of  programs  ranging  from 
the  Saturn  1  and  lb  to  the  Space  Shuttle.  It  has  the  distinct  advantage  in  this 
application  in  that  detail  can  be  applied  where  necessary  and  still  retain  the 
advantage  of  being  able  to  make  modifications  and  changes  in  a  minimum  of  effort. 
Since  the  optical  truss  is  the  most  critical  strqcture,  it  was  modeled  in  detail 
as  a  truss  with  rigid  optical  elements.  The  solar  arrays  were  modeled  as  con¬ 
tinuous  uniform  beams.  The  equipment  section  length  was  assumed  (in  that  no 
data  was  given)  and  modeled  as  a  uniform  beam  with  equal  properties  in  each 
direction. 

It  was  required  to  develop  two  different  models  with  the  same  general  character¬ 
istics  in  order  to  illuminate  features  of  the  control  system  design.  These  two 
models  were  termed  STRENGTH  and  STIFFNESS.  Their  definitions  are  as  follows: 

STRENGTH  MODEL 

The  structural  characteristics  of  this  model  were  envisioned  to  be  representative 
of  a  class  of  vehicles  that  would  be  designed  on  the  basis  of  being  as  light  as 
possible  and  just  being  strong  enough  to  withstand  the  launch  environment.  Thus 
the  strength  model  should  present  low  frequency  modal  data. 

STIFFNESS  MODEL 

The  structural  characteristics  of  this  model  were  envisioned  to  be  as  stiff  as 
possible.  Therefore,  high  frequencies  are  to  be  expected. 

It  is  obvious,  certainly,  that  these  two  definitions  are  not  only  not  very  pre¬ 
cise,  but  are  also  rather  arbitrary.  Therefore,  it  is  very  important  that  the 
developed  structural  model  not  only  meet  the  spirit  of  these  requirements,  but 
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also  be  readily  adaptable  to  fine  tuning  and  adjustments  as 
required.  The  Control  Dynamics  Model  meets  both  of  these  requirements. 

The  various  elements  of  the  Control  Dynamics  model  are  contained  in  the  following 
sections. 
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3.1 


OPTICAL  TRUSS 


The  optical  structure  is  designed  to  be  both  stiff  and  accurate  in  order  to 
maintain  the  critical  dimensions  and  spacing  of  the  optical  elements.  The  study 
model  consists  of  35  nodes  and  117  members.  The  optical  elements  contain  the 
most  significant  mass  of  the  structure  and  are  modeled  as  rigid  bodies  suspended 
on  special  kinematic  mounts  attached  to  the  structure  in  such  a  way  that  no 
structural  loads  are  carried  through  the  fhirrors.  More  will  be  said  on  this 
later. 

There  are  4  elements  in  the  optical  path.  These  are  called  the  primary,  sec¬ 
ondary,  tertiary  and  focal  plane  assembly  respectively.  Since  all  the  mass  is 
concentrated  in  the  optical  elements  which  are  modeled  as  rigid  bodies,  this 
system  has  24  degrees  of  freedom.  Since  there  are  35  nodes  for  which  105  coor¬ 
dinates  are  required  to  specify  the  configuration,  there  are  81  coordinates 
that  must  be  eliminated  from  the  model  description.  Because  of  the  rigor  re¬ 
quired  in  this  elimination,  it  was  considered  important  to  include  a  detailed 
derivation. 

The  positions  of  the  nodes  are  known  from  the  structural  definition.  The  rigid 
body  coordinates  can  be  expressed  in  terms  of  the  nodal  coordinates.  To  do  this 
we  must  first  look  at  the  supports  for  the  optical  elements.  They  are  assumed 
supported  at  three  points  identified  as  A,  B,  and  C  for  convenience.  The  support 
at  point  A  is  considered  to  have  translation  stiffness  which  acts  in  the  three 
coordinate  directions.  At  point  B  the  support  is  considered  to  have  stiffness 
in  only  two  directions  and  provides  no  stiffness  along  the  line  connecting  A  and 
B.  At  point  C  the  support  has  stiffness  only  in  the  direction  normal  to  the 
plane  containing  A,  B,  and  C.  Since  any  or  all  of  the  support  points  can  occur 
at  points  other  than  node  points  of  the  structure,  it  may  be  written  as  a  linear 
combination  of  positions  of  adjacent  nodes.  The  positions  of  the  support  points 
can  also  be  expressed  in  terms  of  the  rigid  body  translation  and  rotation  of  the 
optical  elements. 

The  transformation  between  nodal  and  rigid  body  coordinates  can  be  expressed  as 
follows: 


P  X  =  Q  R,  (1) 

where 

P  -  matrix  to  combine  nodal  displacements  as  required  where  attach  points 
are  located  at  points  other  than  a  node. 
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X  -  vector  of  nodal  displacements 

Q  -  matrix  relating  rigid  body  motion  to  nodal  displacements 
R  -  vector  of  rigid  body  motions  of  optical  elements 

Since  many  of  the  nodal  displacements  do  not  directly  figure  in  the  rigid  body 
motion  of  the  optical  elements,  the  size  of  the  structural  analysis  problem  can 
be  significantly  reduced  by  eliminating  these  coordinates.  This  is  most  easily 
accomplished  by  first  arranging  the  massless  nodes  to  be  the  first  (or  topmost) 
elements  in  the  nodal  position  vector  for  the  system.  The  massless  nodes  must 
be  at  all  times  in  equilibrium  so  that  we  can  write  the  force  vector  in  a  par- 


The  vector  Xy  represents  the  displacement  vector  for  the  massless  nodes.  In  our 
model  78  nodal  displacement  components  do  not  directly  affect  the  position  of 
the  optical  elements.  This  leaves  27  that  do.  Since  the  4  rigid  bodies  require 
only  24  components  to  completely  specify  their  positions,  we  can  perform  a  second 
reduction  step  by  eliminating  these.  The  first  reduction  step  requires  calcula¬ 
tion  of  the  inverse  of  the  upper  part  of  the  structural  stiffness  matrix  K^. 

Since  this  is  a  78  x  78  matrix,  this  step  promises  to  be  an  expensive  step.  This 
process  can  be  simplified  by  taking  advantage  of  the  symmetry  and  spareness  of 
the  matrix  K^.  Since  is  part  of  a  structural  stiffness  matrix  it  is  non¬ 
negative  definite.  This  means  that  the  diagonal  elements  are  nonnegative  and 
are  in  fact  zero  only  if  the  node  is  completely  decoupled  from  the  rest  of  the 
structure  in  that  dimension.  That  is,  if  a  diagonal  element  vanishes  the  whole 
row  and  column  in  which  it  lies  will  also  vanish.  It  is  relatively  easy  to  con¬ 
struct  a  lower  triangular  matrix  Q  which  converts  Ky  into  a  lower  triangular 
matrix  such  that  Q  is  the  identity  matrix.  If  we  define  a  new  vector 

such  that 

Xu  =  Q'^  Xu  (3) 

We  can  solve  for  XJ)  and  eliminate  it  from  the  force  equation  and  determine  an 
effective  stiffness  matrix. 

-  kJ  Q  (4) 

With  this  effective  stiffness  matrix  we  now  wish  to  reduce  the  size  of  our  system 


to  the  minimutn  and  convert  from  nodal  displacements  to  riqid  body  motions  in¬ 
volving  translations  and  rotations  of  optical  elements.  The  reduction  process 
is  similar  to  the  above  in  that  we  select  a  subset  of  the  27  coordinates  to 
eliminate  and  call  them  We  must  solve  for  the  Xy  in  l,erms  of  the  remaining 
coordinates  and  substitute  for  them  in  the  equations  (1)  and  (2)  to  arrive  at  a 
transfomation  expressing  the  remaining  nodal  variables  in  terms  of  rigid  body 
coordinates.  This  transformation  must  then  be  applied  to  the  reduced,  effective 
stiffness  matrix.  The  mathematics  involved  in  this  process  is  straightforward, 
although  tedious,  and  will  not  be  presented  in  detail  in  this  report.  The  re¬ 
sults  of  this  last  step,  though,  is  to  produce  and  effective  stiffness  matrix 
in  terms  of  rigid  body  motions  of  the  four  optical  elements  and  the  stiffness 
properties  of  the  truss  structure. 

Before  presenting  results  for  the  optical  structure,  one  further  observation 
must  be  made.  It  is  noted  that  the  focal  plane  is  mounted  in  such  a  manner  as 
to  allow  essentially  free  motion  along  the  line  of  sight  of  the  optical  truss. 
This  is  considered  an  unsatisfactory  design  feature  and  therefore  was  selected 
as  a  candidate  for  added  stiffness  for  the  stiffness  model.  We,  therefore, 
selected  to  differentiate  between  our  Stiffness  model  and  our  Strength  model 
by  the  stiffness  of  this  section.  The  truss  itself  is  efficiently  designed  to 
carry  the  loads  and  no  appreciable  weight  savings  are  exoected  to  accrue  through 
reductions  in  the  truss  elements.  The  modal  data  for  the  optical  truss  section 
is  presented  in  Table  1  for  the  Stiffness  model  and  in  Table  2  for  the  Strength 
model . 
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3.2 


LIGHT  PATH  ANALYSIS 


The  light  path  analysis  was  accomplished  by  considering  the  mirrors  as  thin 
spherical  mirrors  with  all  centers  of  curvature  on  the  optical  axis.  The 
light  ray  configuration  is  given  in  Figure  2  and  is  based  on  the  referenced 
CSDL  report. 


The  geometry  considerations  utilized  in  determining  the  effects  of  flexibility 
motion  and  rigid  body  motion  at  each  of  the  optical  elements  on  the  light  path 
is  given  in  Figure  3.  We  analyzed  each  mirror  under  the  assumptions  as  given 
and  as  illustrated  in  the  Figure.  The  vector  equation  was  processed  for  each 
mirror  taking  the  object  of  each  mirror  as  the  image  of  the  preceeding  mirror 
in  the  chain  (except  for  the  first  mirror  where  the  object  was  assumed  to  be  at 
infinity).  The  result  of  this  analysis  was  programmed  so  that  we  could  easily 
develop  the  effect  of  each  mode  shape  on  the  light  path  error  (or  apparent  error). 
The  result  is  that  for  each  mode  shape  there  are  three  values  (x,  y,  z  transla¬ 
tion  at  the  focal  plan)  that  represent  the  contribution  of  that  particular  mode 
to  the  error  in  location  of  the  image.  These  three  error  numbers  are  presented 
with  each  mode  in  the  appropriate  section  of  this  report. 
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LIGHT  PATH  ANALYSIS 


Figure  3 


3.3  EQUIPMENT  SECTION 

The  equipment  section  is  composed  of  solar  arrays  and  equipment  module  as  given 
in  Figure  4. 


l.ie  solar  arrays  are  treated  as  uniform  beams  comnencing  at  2  meters  from  the 
centerline  and  extending  to  26  m  on  either  side  in  the  +  X  directions.  The  mass 
of  each  array  is  360  kg.  The  material  and  section  properties  are  as  stated  in 
the  OARPA  reports  with  E  »  1.24  x  lO^^N/M^  ,  A  *  9.407  x  10"^  M^,  I  =  1.874  x 
10*^  M^.  The  central  module  of  the  equipment  section  is  modeled  as  a  uniform 
beam  extending  4m  along  the  negative  Z  direction  starting  at  the  origin  of  coor¬ 
dinates  as  shown  in  Figure  4.  It  has  a  mass  of  3500  kg.  The  rigid  body  proper¬ 
ties  of  the  equipment  section  are  a  mass  of  4220  kg  and  moments  of  inertial  of 
Ixx  =  5890  kg  m2,  =  188069,  I^j  *  185680.  The  three  uniform  beams  are  com¬ 

bined  through  the  process  of  modal  synthesis  to  form  the  equipment  section  model. 
This  resul ts  in  mass  and  stiffness  matrices  as  shown  below 
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The  nodes  that  make  up  the  input  to  the  model  are  developed  as  so-called 
cantilever  modes  having  one  end  clamped  and  one  end  free.  Notice  that  all 
coupling  is  through  the  mass  matrix.  Each  of  the  three  cantilever  beams  of  the 
model  is  truncated  after  the  first  2  bending  modes  in  each  of  the  two  transverse 
axes.  Thus  the  equipment  section  model  has  a  total  of  18  degrees  of  freedom. 

The  material  properties  of  the  center  module  are  El  =  1.788  KIO^  N/Mp  which 
yields  a  first  bending  frequency  of  .5  Hz.  The  solar  array  first  bending 
frequency  is  .38  Hz.  This  model  as  defined  constitutes  our  Strength  Model.  To 
generate  the  Stiffness  Model  the  structures  are  stiffened  by  a  factor  of  10  with¬ 
out  changine  the  mass.  The  modes  of  the  equipment  section  are  shown  in  Tables  3 
and  4. 
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3.4 


COUPLED  ANALYSIS 


The  development  of  the  coupled  model  at  this  point  involves  the  generation  of 
coupled  mass  and  stiffness  matrices  form  the  developed  components.  The  stiffness 
matrix  takes  the  form  as  given  below: 


COMBINED  SYSTEM  STIEPNtSt: 


3^  Oi 
0  2  • 
.  ! 


- br — r 

i 


♦  ?  ICI 

I 


eigenvalues  {natural  frequency)  of  the  optical  truss 
eigenvalues  of  the  equipment  section 
eigenvector  of  the  optical  truss 
eigenvector  of  the  equipment  section 

It  is  noticed  that  the  effect  of  stiffness  changes  in  the  equipment  section 
readily  implemented  in  that  one  would  only  have  to  change  the  value  of  the  eigen^ 
values,  since  the  eigenvector  is  not  changed.  In  addition,  the  coupling  between 
the  equipment  section  and  the  optical  truss  is  apparent  in  the  coupling  terms  in 
the  stiffness  matrix.  Since  the  modal  synthesis  in  this  case  is  acting  on 
uncoupled  free-free  modes,  there  is  no  coupling  in  the  mass  matrix. 

We  therefore  have  all  the  information  that  is  required  to  develop  the  detailed 
coupled  model.  We  have  supplied  the  intermediate  steps  to  facilitate  others 
who  might  like  to  couple  these  two  sections  together  differently.  The  coupling 
was  formed  using  the  form  of  the  isolator  section  as  given  in  the  referenced 
model  by  CSDL.  That  is,  the  isolator  springs  were  considered  to  be  attached  at 
three  points  on  the  equipment  section  and  the  optical  structure. 

The  result  of  this  coupling  operation  is  presented  as  modal  data  in  Figures  5 
and  6.  The  modal  data  is  presented  in  Table  5  for  the  Stiffness  model  and  in 
Table  6  for  the  Strength  model. 


where 
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4.0 


SUMMARY/ CONCLUSIONS 


We  have  presented  a  detailed  structural  model  of  two  forms  of  future  LSS.  The 
form  of  development  of  both  of  these  structures  enables  rapid  and  efficient 
alterations  of  these  models  in  the  future  if  conditions  and  requirements  warrant. 
It  is  thouqht  prudent  at  this  point  to  present  these  models  for  review  and  com¬ 
ment  by  the  VCOSS  contractors  and  others.  The  data  can  be  quickly  changed  to 
emphasize  characteristics  of  LSS  that  are  not  sufficiently  covered  in  these 
current  models. 


A- 14 


TA6LE  1 


PTIFFNESS  MODEL  BENDING  MODES  OF  OPTICAL  STRUCTURE 
(sere  freetuency  inodes  oniitted). 

All  di SF- 1  Bcenients  are  in  nieters. 

Modes  are  nornialised  to  unity  aeneraiized  ness. 

The  followina  table  displays  niodal  diselacenients  •  The  foriiiet 
is  as  shown t 


JOINT  4  JOINT  3  JOINT  6  LOS  ERROR- 


MODE  N  FREQ  =  X.XXXX  HZ 


X  Translation 
Y  Translation 
Z  Translation 

>; «  >:  >t  >; 

V  V*  V* 

#  •  •  *s  y% 

.  V*  \<f 

t\  *\  T% 

M  A  V'  V'  K*  V# 

V#  .  V/  \*  %y  \t 

.  V#  M  M  V 

.  »./  *  /  *.»  ».»  *.» 

/*  •  y\  ts  /%  t  \  f  *  *  \ 

V*  »  V'  yj  M  M  V 

^  *  4  A  #  V  y*  /  *  #  % 

y, .  5-r 

y, .  :  i  t 

>i .  t  ti  iiOii  ty. 

MODE  =  1 

FREQ  = 

2.2092  HZ 

-0.005477 

0.000094 

-0.000321 

0.002359 

-0.003218 

0.015413 

0.002566 

0.003589 

-0.016527 

-0.127966 

-0.005460 

-0.004410 

MODE  =  2 

FREQ  = 

3.5895  HZ 

-0,025537 

-0.000570 

0.000023 

0.051599 

-0.033888 

0.000956 

0.051832 

0.033106 

-0,000788 

0.073843 

0.006635 

-0.001663 

MODE  =  3 

FREQ  = 

5.6413  HZ 

0.001951 

-0,003154 

0.005143 

-0.005546 

-0.003510 

0.005797 

-0,003902 

-0.001794 

0.005326 

-0.004677 

0.190868 

0.009491 

MOPE  ^  4 

FREQ  = 

9,7876  HZ 

0.008534 

0.000664 

0.002317 

-0.012636 

0.009080 

-0.001327 

-0.012658 

-0.007732 

0.002635 

-0.033210 

-0.013211 

0.000962 

MOPE  --  5 

FREQ  = 

10.6198  HZ 

0.001441 

-0,009365 

-0.012824 

-0.002024 

-0.007086 

-0.001294 

-0.002388 

-0,011950 

-0,000018 

-0.006785 

0,123783 

0.004621 
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MODE 

= 

6 

FREO 

S 

-0.021381 

0 

0.001019 

-0 

-0.002751 

0 

MODE 

= 

7 

FREO 

= 

-0.001181 

0 

-0.000207 

-0 

-0.004042 

-0 

MODE 

:s 

8 

FREO 

= 

0.000851 

-0 

0.002901 

0 

0.019083 

0 

MODE 

9 

FREO 

= 

-0.00D240 

-0 

0.001101 

0 

0.001115 

-0 

MODE 

10 

FREQ 

s 

-0.000569 

-0 

0.003696 

0 

0.008617 

-0 

MODE 

11 

FREO 

= 

0.026150 

-0 

0.009632 

0 

-0.031006 

-0 

MODE 

12 

FREO 

= 

-0.005105 

-0 

0.004906 

0 

0.004028 

-0 

11.6606  HZ 


038098 

024827 

007606 

0.038327 

0.027090 

-0.008909 

-0,006099 

0.000001 

-0.004296 

18.3524 

HZ 

005581 

002376 

001821 

0.005063 

0.001783 

0.003761 

-0.011199 

0.015729 

0.004589 

19.0918 

HZ 

010937 

006764 

003554 

-0.009135 

-0.000437 

-0.000746 

0.039565 

0.040702 

-0.025179 

23.7055 

HZ 

004564 

003056 

000186 

-0.004179 

-0.000954 

-0.004018 

0.019475 

0.072887 

-0.032830 

24.4272 

HZ 

007718 

006760 

004658 

-0.007326 

0.000474 

-0.010399 

0.034363 

0.063281 

-0.039389 

30.5696 

HZ 

057293 

051865 

005163 

-0.059287 

-0.034757 

-0.005794 

0.051965 

-0.056568 

0.003104 

32.5465 

HZ 

005037 

004464 

007363 

-0.004503 

0.004884 

-0.016295 

0.037886 

-0.297196 

0,035895 
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MO  TIE  =  13 


EREO  = 


36.0537  H2 


-0.004611 

0.013281 

0.012046 

0.076989 

-0.005616 

-0.009569 

-0.000164 

-0.008783 

0.016966 

0.011286 

-0.007370 

0.005561 

b 


MODE  *  14 

FREQ  * 

52.3519 

HZ 

0.042057 

0.028234 

0.025860 

-0.633964 

0.008891 

0.017854 

■0.001896 

-0.067880 

0.029311 

-0.09536B 

0.071788 

0.039530 

MODE  =  15 

FREQ  = 

66.7940 

HZ 

0.015749 

-0.017362 

0.019001 

0.604286 

‘ii 

0.000231 

-0.013148 

0.016127 

-0.232495 

0.038447 

0.040765 

0.128205 

0.041543 

r  • 

1  • 

i-'-- 

MODE  =  16 

FREQ  = 

71.6937 

HZ 

0.020971 

-0.039879 

-0.019135 

-0.002674 

-0.046937 

0.148159 

-0.002577 

-0.031331 

0.117787 

0.107122 

0.704694 

-0.008513 

MODE  =  17 

FREQ  * 

83.0792 

HZ 

0,039134 

0.025936 

0.006834 

-0.016030 

0.034106 

0.045947 

-0.012028 

0.017050 

-0.001235 

0.218509 

0.130179 

0.015092 

MODE  -  18 

FREQ  = 

96.3105 

HZ 

0.036327 

0,016169 

-0.062904 

0.030580 

0.007441 

-0.074183 

0.040541 

0.030574 

-0.064852 

-0.008927 

-0.393105 

-0.057276 

r. 


Opticel  Structure  Rigid  Body 


Properties  * 


Mass  -  5000.0  Ks 


Center  of  mess 


0  t  1.5  f 


14  )  n. 


Inertia  Matri;: 

640500 

0 

0 

about  Mass  Center  = 

0 

530000 

-10000 

ka(ni»*2) 

0 

-10000 

210500 

h 
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TABLE  2 


STRENGTH  MODEL  BENDING  MODES  OF  OPTICAL  STRUCTURE 
(zero  freouency  modes  omitted). 

All  disF-lacements  ere  in  meters. 

Modes  ere  normelized  to  unity  Senerelized  mess. 

The  following  table  displays  modal  displacements.  The  format 
is  as  shown. 


JOINT  4 

JOINT  3  JOINT  6 

LOS  ERROR 

MODE  =  N 

FREQ  = 

X.XXXXZ 

X  Translation 

Y  Translation 

Z  Translation 

X  «  X  X  )■%  X 

A  4  A  A  A  A  A  •  V  4  « 

A  4  A  A  »s  A  A  A  4  A  »  »  A  A  A  A 

X  4  X  X  J :  )■: 

X  4  X  X  J-C  X  X 

X  4  X  :  X  X  x  x 

MODE  =  1 

FREQ  = 

0.2494HZ 

0.000003 

0.001009 

-0.00(£.684 

0.000022  0.000019 

0.001008  0.001008 

-0.007402  -0.007391 

-0.000018 

-0.193395 

-0.004965 

MODE  =  2 

FREQ  = 

2,1996HZ 

0.005536 

0.000102 

0.000165 

-0.001941  -0.002220 

0.003290  -0.003284 

-0.015536  0.016395 

0.128133 

0.002569 

0.004130 

MODE  =  3 

FREQ  = 

3.5895HZ 

-0.025527 

-0.000414 

-0.000244 

0.051836  0.052012 

-0.033762  0.033270 

0.000712  -0.001091 

0.073652 

-0.000852 

-0.001965 

MODE  =  A 

FREQ  = 

9.7872KZ 

O.OOB470 

0.000716 

0.002384 

-0.0J2603  -0.012634 

0.009140  -0.007698 

-0.001325  0.002630 

-0.033130 

-0.014812 

0.000885 
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MOHE  = 

5 

FR-EQ  = 

10.5 

443HZ 

p 

0.002368 
-0 . 009401 
-0.012620 

-0.002883 

-0.007190 

-0.001341 

-0.003009 
-0.011705 
0.00011 1 

-0.007381 

0.134121 

0.005554 

MODE  « 

6 

FREQ  = 

11.8486H7 

-0.021352 

0.001101 

-0.003063 

0.038412 

-0.024471 

0.007573 

0.038439 

0.026781 

-0.008876 

-0.006034 

-0.000846 

-0.004644 

MOPE  ^ 

7 

FREQ  = 

ie.3l76H2 

1 

0.002021 

0.002119 

0.006853 

-0.005283 

0.004958 

0.002292 

-0.005126 

-0.000702 

-0.003791 

0.017235 

-0.010512 

-0.009683 

MOPE  = 

8 

FREQ  = 

-0.001950 

-0.005643 

-0.016733 

18.75 

0.005435 

-0.009650 

-0.002980 

60H2 

0.005334 

-0.001450 

-0.000298 

-0.035037 
-0 . 038066 
0.025622 

MOPE  * 

9 

FREQ  * 

23.6388H2 

* 

0.000496 

0.003502 

0.001082 

-0.000909 

0.004307 

-0.001315 

-0;001249 

0.002389 

-0.005282 

0.018886 

0.072333 

-0.036225 

! 

MOPE  = 

10 

FREQ  = 

24 . 3000H2 

■'» 

0.000403 

0 . 005935 
0.006996 

-0.001191 

0.006849 

-0.005202 

-0.001862 

0,004601 

-0.0095B8 

0.025116 

0.042542 

-0.032689 

MOPE  = 

1 1 

FREQ  = 

28 . 4961  HZ 

>4 

-0.010490 

0.005302 

0.011926 

0.051218 

-0.028396 

-0.003073 

0.050314 

0.039922 

0.000668 

-0.056870 

-0.067494 

0.013315 

r 


A- 19 


Of^ticel  Structure  Rigid  Poda 


Properties  J 

Mess  =  5000.0  Kg 

C e ri t e r  of  ni e s s  =  (  0  »  1.5  »  14  )  iti 

InertieMatri:;  640500  0  0 

about  Maso  Center  =  0  530000  -10000  Kg(ni*f2) 

0  -10000  210500 


TA&LE  3 


STIFFNESS  MODEL  DENDING  MODES  OF  EOUIFMENT  SECTION 
All  disFlocenier. ts  sre  iri  nieters. 

Modes  ere  riornielized  to  uriity  serierslired  i\ts 

The  following  table  disf-lavs  nodal  di  s^  1  acenente  .  The  fo 
is  as  shouri. 


JOINT  42  JOINT  4^-  JOINT  47 

MODE  =  N  FREQ  =  X.XXXX  HZ 

X  Translation  ; ;  : ; . : ;  :<  :<  :<  :■ 

Y  Translation  •  /"j  <■*  j*?  ^  m  #  x  **»  !■»  y*  * 

1  T  r  a  n  s  1  a  t  i  0  ri  :  : : ;  : ;  ;  : :  :<  .  ;< :  :  :< ; 

X  Rotation 
Y  Rotation 
2  Rotation 


ACTUATOR 


MODE  = 


FREQ  = 


0.0000  HZ 


0 .01539373 
0.00000000 
0 . 00000000 


0.01539373 
0.00000000 
0 . 00000000 


0.0153?3“3 
0. 00000000 
0.00000000 


0.01539373 
0.00000000 
0 . 00000000 
0 . 00000000 
0 . 00000000 
0.00000000 


MODE 


FREQ  = 


0 . 00000000 
0.0153®373 
0 . 00000000 


0.00000000 

0.01539373 

0.00000000 


0.0000  H7 

0. 00000000 
0. 015393“3 
0.00000000 


0 . 00000000 
0.01539373 
0 . 00000000 
0 . 00000000 
0 . 00000000 
0 . 00000000 


MODE 


FREQ  = 


0.0000  HZ 


-0 . 00000000 
-0 . 00000000 
0. 01 53«373 


-0 . 00000000 
-0.00000000 
0.015393~3 


-0 . 00000000 
-0.00000000 
0.01539373 


-0 . 00000000 
-0. 00000000 
0 . 0 1  5  3  3  7  3 
0 . 00000000 
0 . 00000000 
0.00000000 


MODE  =  4 


FREQ  =  0.0000  HZ 


0 00 00000 
0 , 0  :  1  e  1  6  6  7 
0 . 0  ^  5  1  4  9  6  9 


0 . 00000000 
-0 . 021 6 1 6^  ■ 
-0. 0e51  4  =  6'=' 


0.0000  .'00 
-0.0216:66? 
-0.06514=^9 


0 . 00000000 
-0.01015032 
0 . 0  c-  0  C  0  0  0  0 
0.013C2®‘=4 
0 . 00000000 
0 . 00000000 


c  e  , 

mat 


SENSOR 


0.01539373 
0 . 00000000 
0 . 00000000 
0 . 00000000 
0 . 00000000 
0.00000000 


0 . 00000000 
0.01539373 
0 . 00000000 
0 . 00000000 
0.00000000 
0. 00000000 


-  V  .  V'  V 0  V  V 0 
-0.00000000 
0 . 0153=373 
0 . 00000000 
0 . 00000000 
0. 00000000 


0 . 00000000 
0 . 0305030P 
0 . 00000000 
0 . 01302=94 
0 . 00000000 
0.00000000 
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.  FREQ 


0.0000  HZ 


0 . 00382550 
•0. 00000000 
0 . 00000000 


■0.01160346 

•0.00000000 

0.00000000 


■0.00000003 

0.00000004 

0.00137836 


0.00000005 

0.00412883 

0.00825548 


0.00096677 
0.00000006 
0. 00000014 


0. 002?4®e.3 
0.00000009 
0.00000018 


0.00382550 
-0 . 00000000 
0 . 00922362 


0.011 60346 
-0.00928277 
-0.00000000 


-0 . 00000002 
0 . 00000003 

0 . 00137849 


-o.oooooooe 

-0.00412378 

-0.00825549 


0 . 00096628 
0 . 00000014 
0 . 00232868 


-0 . 00294970 

0 . 00235^86 
-  0 . ;  000001 


0.00382550 
-0 . OOOCOOOC 
-0 . 0092236  2 


0.0000  HZ 

0 .01160346 
0.00928277 
-0.00000000 


1.2140  HZ 

-0 . 00000002 
0 . 00000004 
0.001  3785=' 


1 .2238  HI 

-0.00000008 
-0 . 00*1 1  2388 
-o.ooe:'5542 


1.2472  HZ 

0 . 00096628 
-0 . 00000026 
-0 . 00232885 


1 .2477  HZ 

-0.00294979 
-0 . 00235968 
0. 00000038 


0.00179630 
-0 . 00000000 
-0  .  OOOO'OOOO 
0 . 00000000 
0.00230590 
0.00000000 


0.00000000 
-0 . 00000000 
0.00000000 
0 . 00000000 
0 . 00000000 
0.00232069 


0.00000042 
0.00000101 
0.00137845 
-0 . 00000440 
0 . 00000264 
0 . 00000000 


0.00000023 
-0.00261930 
0.00000001 
0.00131489 
-0.00000012 
-0 . 00000001 


0. 00032665 
0. 00000015 
0 . 00000003 
-0 . 00000081 
-0. 00001432 
-0 . 00000005 


0.00  0  0  0  0  0  ~ 
0. 00000008 
0 .00000000 
-  0 . 0  0  0  0  0  0  0  3 
-0 . 00000004 
-0. 00  0  5899-3 


FR'EO  = 


FREQ  = 


FREQ  = 


FREQ  = 


FREQ  = 


-0.00539812 
0 . 00000000 
-0 . 00000000 
0 . 00000000 
0 . 00230590 
0 . OOOOOOOO 


-0 . OOOOOOOO 
0.00000000 
0.00000000 
0.00000000 
0 . OOOOOOOO 
0 . 00232069 


0,00000745 
-0 . 00001273 
0 . 00137845 
-0 . 00000925 
0 . 00000489 
0 . OOOOOOOO 


-0 . 00000006 
0.00160497 
0.00000001 
0.00102253 
0 . 00000095 
-0.00000001 


-0 . 00306323 
-0. 00000240 
0 . 00000003 
-0 . 00000177 
-0.0  C'  ■:■  6  5582 
-0 . 00000005 


0. 0  0000025 
-0 . 0000001  5 
0  .  OOOOOOOO 
-0  .  000000 1  2 
0 . 0  0  C'  0  0  0  4  7 
-0.00058994 


MOt'E 


FftEQ 


1.5817  HZ 


1  1 


0.00037809  0.00037808 

-0.00000262  -0.00000262 
0.00003619  -0.00008029 


0.00037808 
0 . 00000262 
0 . 00000795 


-0.0761 1863 
-0.0018172® 
0.00000001 
0 . 00c  3  ■120® 
-0 . 35018975 
-0.00000000 


-0 .99964207 
0 . 02382944 
0 . 0000000  1 
0.01748096 
-0.7341 3873 
-0 . 00000000 


MODE  =  12 


0 . 00000899 
0 . 0001 1005 
-0.00152032 


FREO  = 

0 . 00000899 
0 . 0001 1005 
0.00151923 


1.5821  HZ 

0 . 00000899 
0.00011005 
0.00152134 


-0.00181123 
0 .07637052 
-0 . 00000002 
-0.35057145 
-0.00833309 
0 . 00000000 


-0 . 02373733 
-1 .00141907 
-0 . 00000002 
-0 . 73462784 
-0.01746915 
0 . 00000000 


MODE  •-  13 


mode  =  14 


M  2  .'I  E  -  15 


mode  =  l9 


0 .00000190 
-0 . 00000022 
0.00071949 


0.00020668 

0.00000001 

0.00000261 


-0 . 00060411 
0 . 00000006 
0. 00000304 


0 . 00000002 
0 . 0021  76  26 
-0 . 00435272 


FREQ  = 

-0 . 00000338 
0 . 000001S9 
0 . 00071681 


FREQ  * 

0.00019031 

0.00000656 

0.00048152 


FREQ  - 

0.00060946 
-0 . 00048537 
0 . 0000097-5 


FREQ  = 

-0 .C0000002 
0. 00217628 
0 . 00435287 


7.5861  HZ 

-0.00000338 
-0 . 00000234 
0.00072038 


7.5875  HZ 

0.00019031 
-0.00000654 
-0 . 00047620 


7.5876  HZ 

C  .  00060946 
0.00048549 
-0 . 0000031 6 


9, *028  HZ 

-0 . 00000002 
0.00217625 
0.00435287 


-0.00000133 
-0 . 00000012 
0 . 00071904 
-0 . 00000027 
0 . 00000052 
-0.00000053 


0.00015214 
0. 00000001 
0 . 00000264 
0 . 00000001 
0 . 0001372® 

-0 . 00000164 


0 . 00000201 
C . 00000009 
0 . 0  0  0  0  0  3  1  ® 
0 . COOCCOO' 
0 . 00000161 
0 . 00012136 


0.00  C  0  0  0  ® 1 
0 . 0  0  2  C'  6. 1  2  6 
0 . 0  0  0  0  0  C'  0  7 

-0 .  cooc9;'®o 

-0.00000143 

-  0 . 0  0  0  0  0  0  0  c 


0 .00000270 
-0 . 00000059 
0.00071904 
-0 . 00000041 
-0 . 000003S1 
-0.00000053 


-0 . 00017924 
0.00000001 
0 . 00000264 
0 . 00000003 
0 . 00040573 
-0.00000164 


-0 . 00000309 
-0.00000015 
0 . 00000317 
-0  .  OOO'OOOI  5 
0 . 00000500 
0.00012136 


-0 . 00000295 
-0.00344018 
0 . 00000007 
-0 . 002 9 48 75 
■  0 . 0  0  0  0  0  2  5  1 
-0.00000000 


A-24 


FREQ 


9.9097  HZ 


1  7 


-0.00021255 
-0.00000012 
0 . 00000018 


-0 . 00021255 
-0 . 00000012 
-0 . 00001642 


-0.00021255 
-0 . 00000012 
0.0000160'’ 


0 . 347568®6 
0 . 00009289 
0.00000000 
0.00013318 
-0.49762329 
0 . 00000000 


-0 . 99994171 
-0.00026724 
0.00000000 
-0.00027223 
1 . 01880682 
0 . 00000000 


18 


FREQ  =  9.9102  HZ 


-0.00000006 

0.00042336 

-0.00063901 


-0 . 00000006 
0.00042336 
0.00063900 


-0.00000006 

0.00042336 

0.00063901 


0.00009293 

-0.34746578 

0.00000000 

-0.49784425 

-0.00013306 

-0.00000000 


-0 . 00026738' 
1 . 00017655 
0 . 00000000 
1 .01856959 
0 . 00027240 
-0.00000000 


EouiFnier. t  Sectior.  Risiii  Body 
F'rop'ertie;  3 

Hass  =  4220.0  Fs 

Center  of  ni3;s  =  <  0  »  0  t  -1.659  >  n 

Inertia  Matri;;  5&90  0  0 

bout  Hass  Center  =  0  188069  v  3'S(ti.»*2'' 

0  0  185680 


TABLE  4 


X 

- 

J| 

■, 

- 

X 


JTRE^^GTH  model  B'ENDING  modes  of  EOL'IRMENT  section 

All  displace rtiente  are  ir,  n'etef's^ 

Modes  are  ri or fii sliced  to  uriitv  serie'aliced  niass. 


The  followiris  table  display's  n.odal 
IS  as  -showri. 

d  1  s  r- 1  a  c  e  Hi  e  ri  t 

s .  The  To 

r  fii  a  t 

JOINT  42  JOINT  46 

JOINT  47 

ACTUATOR 

SENSOR 

MODE  =  N 

FREQ  = 

X.XXXX  H2 

X  Trarislatiori 
Y  Trarislation 
Z  Trarislatior; 

«•»  4  »•*  I'J  *  »  »  •  V  '  w  »  .  »  »  t  /  %  »•» 

•  J  !  j  J  i  t  i  t  i  t  X  i  t  i  t  •  i  t  i  t  i  t  i  t  i  t  i-t  i  t 

it  *  it  it  i  t  i't  it  i  >  i't  i'»  i\  •  t  %  *\ 

i*t  •  X  i't  i't  i  t  #•.  /  « 

i'.  ♦  <  t  »■»  «'♦  t  .  ./•»<» 

X  Rotation 

Y  R  0 1  a  1 1  0  ri 

Z  Rotatiori 

it  •  i  t  i  t  i  t  >  \  i'»  /  <  r\  t\ 

t'\  •  /•  i  t  /•*  /•*  i  t  f  % 

i  t  «  i  t  i  t  i  t  i  t  i  t  i  t  i  t  i't 

it  •  i't  it  /'t  i'«  >  k  i'«  /  >  /  t 

it  •  it  i  t  / 1  i  t  >  •  i  t  i  t  / « 

i't  *  i  t  i  «  i  t  i  t  i  t  i  t  i  t  i't 

MODE  =  1 

FREQ  = 

0.0000  HZ 

0 .01539373  0 .01539373 
0.00000000  0.00000000 
0.00000000  0.00000000 

0 .01539373 
0.00000000 
0.00000000 

0 . 01539373 
0 . 00000000 
0 . 00000000 
0 . 00000000 
0 . 00000000 
0 . 00000000 

0 . 01539373 
0 . 00000000 
0 .00000000 
0 . 00000000 
0 . 00000000 
0 .00000000 

MODE  ^  2 

FREQ  = 

0.0000  HZ 

0 . 00000000  0 . 00000000 

0 .01539373  0. 01 539373 

0 . 00000000  0  .  OCOOOOOO 

0.00000000 

0.01539373 

0. 00000000 

0 . 00000000 
0.01539373 

C’ .  0  0  0  0  0  0  0  0 
0 . 00000000 
0 . 0  000  0  00  0 
0 . 00000000 

0 . 0000:  000 
0.01539373 
0 . 0  0000  00  0 
0 . 00000000 
0 . 00000000 
0 . 00000000 

MODE  =  3 

FREQ  - 

0.0000  HZ 

-  0  .  OOC'OCOOO 
-0  .  OOC'OOOCO 
0 . 01539373 


-0 . 00000000 
-0 . 00000000 
0 . 0  1  5  3  ®  3  7  3 


'0 . 00000000 
-0.00000000 
0.01539373 


-  ■:  .00000  0'  0  0 
-C . 00000000 
0 .01539373 
0 . 00000000 
0 . 00000000 
o.ocoooooo 


-0 . 00  00  00  0  0 
-0 . 00000000 
0 . 0  1  5  3  =■  3  7  3 
0 . 00000000 
0 . 00000000 
0 . 00000000 


MODE 


4 


FREO  = 


0.0000  HZ 


0 .00000000 
-0. 021^1667 
0.06514969 


0 . 00000000 
-0 . 02161667 
-0.0651 4969 


0.0000  0  0  0  0 
-0.02161667 
-0.06514969 


0.00000000 
-  0 . 0  1  0  1  5  0  3  2 

0 .000'  V  0  000 

0.01 302994 
0 . 00000000 
0 . 00000000 


0 . oooooooc 

0 . 03050308 
0 . 0  0  0  0  0  00  0 
0.01302994 
0 . 00000000 
0 . 00000000 


I 
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1 


MODE 


fREQ 


0,0000  HZ 


MODE  * 


MODE  *■ 


MODE  = 


MODE  =■ 


MODE  = 


0.00J82550  0.00382550  0,00382550  0.00179630  -0.00539812 

0.00000000  -0.00000000  -0,00000000  -0,00000000  0.00000000 

0.00000000  0.009223-62  -0.0092236-2  -0.00000000  -0.00000000 

0.00000000  0.00000000 

0.00230590  0.00230590 

0.00000000  0,00000000 

FREO  =  0.0000  HZ 


-0 .01 1 60346 
-0 . 00000000 
0.00000000 


0.011 60346 
-0 . 00928277 
-0 . 00000000 


0.011 60346 
0 . 00928277 
-0.00000000 


0.00000000 

-0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00232069 


-0 . 00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00232069 


7 


FREO 


0.3839  HZ 


-0.00000012 
0. 00000023 
0.00137800 


0.00000003 
0 . 0000001  7 
0.00137880 


0.00000003 
0 . 00000029 
0 . 00137901 


0 . 00000008 
-0. 00000037 
0,00137845 
0.00000179 
0 . 00000328 
0 . 00000001 


0.00000893 
0 . 00000584 
0.00137845 
0 . 00000436 
0.0*000412 
0 , OOOOOOOl 


8 


FRED 


0.3870  HZ 


-0 . 00000014 
0 . 0041 2883 
-0.00825554 


0.0000001  1 
0.00 11 2873 
0 .00825536 


0 . 00000011 
0, 0041 2893 
0 . 00825542 


0.00000054 
0.00261919 
-0.00000008 
-0.00131405 
0.00000012 
0 . 00000002 


0.00000085 
-0. 00160293 
-0.00000008 
-0.00102108 
0 . 00000259 
0 . 00000002 


9  FREO==  0.3944HZ 

-0.00096867  -0.00096439  -0.00096439  -0.00032582 
-0.00000005  -0.00000176  0.00000166  -0.00000120 

0.00000006  -0.00232895  0.00232858  -0.00000006 

0.00000516 
0 . 00002038 
0.00000043 

10  FREQ  =  0.3945  HZ 

-0.00  294901  0.0  0  295041  0.00  2'?  5  041  0.00  0  000  6  0 

-0.00000017  -0.00235994  0.0023595®  -0.00000044 

0.00000032  0.00000130  -0.00000208  -0.00000003 

0.00000166 
0 . 0000009S 
0 , 00058994 


0 . 0030B009 
0 . 00001494 
- 0 . 00000006 
0.00001098 
0.00066633 
0 . 00000043 


0 . 00000074 
0 . 00000456 
-0 . 00000003 
0.00000333 
0 . 0  0  0  0  0  2  2  6 
0 , 00058994 


A-27 


MODE  ^ 


MODE  = 


MODE  = 


MODE  = 


MODE  = 


MODE  = 


1  1 


EREG  =  0.5002  HZ 


0.00037815  0.0003'^ei7 

-0.00000088  -0.00000089 
0.00001226  -0.00005641 


0 . 00037817 
-0 . 00000087 
0 . 00003192 


-0.07613716 
-0.00061613 
0 . 00000001 
0 . 00282901 
-0 . 35027820 
0 . 00000000 


-0 . 99989396 
0 .00808035 
0.00000001 
0 . 00592754 
-0.73432106 
0 . 00000000 


12 


FREQ 


0.5003  HZ 


0.00000302 
0 . 00011008 
-0.00152070 


0.00000305 

0.00011007 

0.00152032 


0.00000305 

0.00011009 

0.00152110 


-0.00061372 

0.07638976 

0.00000001 

-0.35065913 

-0.00282648 

0.00000000 


-0.00806775 
-1 . 00167036 
0.00000001 
-0.73481214 
-0.00592246 
0.00000000 


13 


FREQ 


2.3989  HZ 


0.00000063 
-0 . 00000039 
-0 . 00071828 


-0.00000019 

-0.00000005 

-0.00071930 


-0.00000019 

-0.00000072 

-0.00072035 


-0 . 000001  10 
-0 . 00000058 
-0.00071906 
-0 . 00000021 
0 . 00000166 
-0 . 00000008 


0.00000263 
0.00000140 
-0.00071906 
0. 00000127 
-0.00000361 
-0.00000008 


14 


FREQ  =  2.3994  HZ 


0.00020708 
0.00000002 
0 . 00000076 


0.00018991 
0.00000689 
0 . 00047969 


0.00018991 

-0.00000685 

-0.00047802 


0.00015235 
0 . 00000005 
0 . 00000080 
0.00000000 
0 . 00013680 
-0 . OOOCOl 72 


-0.00018037 
-0.00000019 
0 . 00000080 
-0.00000016 
0 .0004061  1 
-0 . 00000172 


FREQ 


2 . 3994  HZ 


-0.00060398  0.00060960 
-0.00000018  -0.00048561 
-  0 . 00000012  0 . 00000594 


0 . 00060960 
0 . 00048525 
-0.00000761 


0.00000171 
-0 . 00000022 
-0 . 00000048 
-0 . 00000004 
0 . 00000213 
0.00012136 


-0 . 00000245 
0 . 00000046 
-0 . 00000048 
0 . 00000041 
0.00000396 
0. 00012136 


16 


FREQ  = 


2.4042  HZ 


0 . 0  00  0  0  0  0  5 
-0 . 0021 7626 
0.00435292 


-0 . 00000005 
-0.00217622 
-0 . 00435267 


-0.00000005 
-0. 0021 7630 
-0 . 00435267 


0 . 00000071 
-0. 002061 15 
0.00000013 
0. 0000930S 
-0.00000110 
-0 . 00000001 


-0 . 00000226 
0. 00343988 
0 . 00000013 
0 . 00294843 
0 . 0000C201 
-0 . 00000001 


A-28 


MODE 


FPEG 


3. 1337  HZ 


1  7 


R 


0.00021255  0.00021255 

-0.0000004S  -O.OOOOOOdS 
0.00000102  0.00001523  -0.00001727 


0.00021255  -0. 34756863  0.9999788*- 

0.00000068  0.00055795  -0.00160716 

0.00000000  0.00000000 
0.C007'?876  -0.0016363: 
0.49762210  -1 .01880670 
O.OOOOOOOO  0.00000000 


MODE  = 


18 


.  00000024 
.  00042336 
i  00063901 


FREQ  = 

0.00000034 
0.00042336 
0.00063903 


3.1339  HZ 

0.00000034 

0.00042336 

0.00063898 


-0.00055803  0.00160501 

-0.34746537  1.00017500 

0.00000000  0.00000000 

-0.49784374  1.01856816 

0.00079879  -0.00163589 

0.00000000  0.00000000 


Ec»<jiPn.er.t  Sectior.  Rigid  Dody 


Prop-ertles 

Mess 

Ceriter  of  n, ess 

Inertia  Metn:; 

?  b  0  'j  t  Mess  C  e  ri  t  e  r 


4220.0  kg 

(  0  .  0  .  -1.659  )  n, 

5890  0  0 

0  188069  0 

0  0  185680 


k  S  ( It. »  »  2  ) 


1 
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table  5 


STIEFNESS  mo  I' el  ben  BING  MODES  OF  COMBINED  SYSTEM  MODEL 
(zero  freoijer.cy  itrodes 
All  d  1  £  P- 1  5  c  e  (»i  e  n  i  s  in  meters. 

Modes  ere  normelized  to  units'  senerelized  mess. 
The  followins  teMe  displays  modal  displacements.  The  format 


is  3S  shown 

ACTUATOR 

SENSOR 

LDS  MOTION 

MODE  =  N 

FREQ  = 

X.XXXX  HZ 

>;  Translation 

Y  Translation 

Z  Translation 

X  Rotation 

Y  Rotation 

Z  Rotation 

»'*  • 

y*  t  t't  *'%  /'«  y* 

t'\  t  < « «'« f's  y* 

y\  *  /«  «*\ «'« *'« 

V#  «/««»«\#  \y  »,/ 

t%  « 

\  4  \  \  \ 

/  •  t  ^  » 

^  *'»  *\  •*•  ^*  a* 

M  .  ».»  V  *•' 

*■!  a  ;<  ;<  j  i  )  ■;  >i  >i 

;; . 

MODE  =  1 

FREQ  = 

0.1772  HZ 

-0 . 003313 
-0.000005 
0.000014 
-0.000000 
-0.00210S 
0 . 000074 

0 . 002764 
-0.000003 
0.000014 
0.000002 
-0.002112 
0.000074 

-0.000001 

0.000294 

-0.001067 

mode  *  2 

FREQ  * 

0.3552  HZ 

-0 . 002630 
-0 .000003 
-0 . 000021 
-0.000002 
-0.000126 
0.001623 

-0.002736 
-0 . 000009 

-o.oooon 

-0.000002 

-0.000218 

0 . 001623 

-0.000010 

0.000005 

0.000013 

mode-  =  3 

FREQ  = 

0.4040  HZ 

-0 . 000007 
0 . 001 085 
0.010924 
0.001403 
0 . 000000 
0 . 000003 

-0.000008 
0 . 005469 
0 . 010924 
0.001226 
-0 . 000001 

0 . 000003 

0 . 000098 
0.00001 1 
-0.000003 

mode  =  4 

FREQ  = 

0.5304  HZ 

-0 . 000027 
0 . 007635 
-0.002575 
-0 . 000221 
-0.000003 
-0 . 000001 

-0.000029 
0.006‘’16 
-0 . 002575 
-0.001100 
-0. 000006 
-0.000001 

-0 . 005480 
0.001432 
0.000747 

A- 30 


MOI'E 


FREQ 


0.6196  HZ 


0 . 008122 
0.000021 
-0.000004 
0. 000000 
0.001025 
0.000435 


0 . 009431 
0.000021 
-0.000004 
-0 . 000004 
0.001941 
0.000435 


-0 . 001432 
0.001054 
0 . 000000 


MODE 


6 


f'REQ  =  1.1925  HZ 


-0.000000 
-0.000898 
-0.000361 
-0.000141 
0.000001 
0 . 000000 


0.000002 

-0.001116 

-0.000361 

-0.003883 

0.000001 

0.000000 


0 . 009431 
0.000021 
-0 . 000004 


MOPE  =  7  FREQ  =  1.2144  HZ 


-0 , 000000 
0.000137 
-0.001487 
-0.000093 
-0 . 000002 
0.000000 


-0.000005 
-0.000154 
-0. 0014S7 
-0.000201 
-0.000003 
0.000000 


-0 . OOOOOP 
-0. 124940 
-0.001087 


MOPE 


8 


-  1.2480  HZ 


-0.000418 
0 . 000000 
-0 . 000000 
0. 000001 
-0 . 000050 
0.000058 


0.002853 

0.000002 

-0.000000 

0.000001 

0.000529 

0.000058 


0.000003 

-0,000085 

0.000018 


MOPE  =  9 

-0 . 000140 
0 . 000000 
-0 . 000000 
0.000000 
-0 . 000082 
-  0 . 000620 


FREQ  =  1.2499  HZ 

-0.000007  -0.000271 

-0.000000  -0.000036 
-0.000000  -0.000001 
0 . 000000 
-0.000107 
-0.000620 


A-31 


MorT  = 

10 

FREQ  = 

1.4901  HZ 

-0.000001 
0 . 000294 
-0.001087 
-0.041142 
-0.000003 
-0.000000 

-0.000008 
-0.124940 
-0.001087 
-0 . 099322 
-0.000007 
-0.000000 

-0.005480 

0.001432 

0.000747 

MODE  = 

1 1 

FREQ  = 

1.5817  HZ 

-0.076121 
-0.000274 
-0.000001 
0.001239 
-0.350286 
0 . 000001 

-0.999901 

0.003533 

-0.000001 

0.0025B7 

-0.734339 

0.000001 

-0.000190 

-0.000297 

-0.000525 

MODE  = 

12 

FREQ  = 

1.5846  HZ 

-0.000270 

0.077243 

0.000172 

-0.348482 

-0.001245 

0.000000 

-0.003552 

-0.994330 

0.000172 

-0.728201 

-0.002609 

0.000000 

0.000000 

0.000000 

0.000000 

MODE  = 

13 

FREO  = 

2.2282  HZ 

-0.000015 
0.000029 
0.000012 
-0.000069 
-0 . 000063 
-•0.000012 

-0.000056 

-0.000182 

0.000012 

-0.000111 

-0.000110 

-0.000012 

0.009431 

0.000021 

-0.000004 

MODE  = 

14 

FREQ  = 

3.7563  HZ 

0 . 000184 
-0 . 000007 
0 . 000001 
0 . 000003 
-0.000144 
0 . 000055 

-0.000245 

0.000006 

0.000001 

0.000005 

-0.000250 

0.000055 

0.000001 

0.000529 

0.000058 

A- 32 


MODE  =  15 


MODE  =  16 


MODE  =  1 7 


-0.000010 
0.000005 
0.000018 
-0.00002 
0 . 00000 
-0.00000 


-0 . 000001 
-0 . 000002 
0.000720 
-0.000001 
-0.000000 
-0.000001 


0.000152 
"0 . 000000 
0.000003 
-0.000000 
0.000137 
-0.000006 


EQ  = 

0.000003 

-0.000085 

0.000018 

-0.000046 

-0.000001 

-0.000001 


EREQ  = 

0 . 000001 
-0.000003 
0.000720 
-0.000001 
-0.000004 
-0.000001 


FREO  = 

-0.000180 

0.000001 

0.000003 

0.000001 

0.000405 

-0.000006 


MODE  = 


MODE  = 


18 


FREQ  « 


0.000008 
0 . 000000 
0 . 000004 
0 . 000002 
0 . 000005 
0 . 000121 


-0.000014 
-0.000001 
0 . 000004 
-0,000001 
0.000019 
0.000121 


FREO  =  7 


0 . 0  0  C  0  0  1 
0.002112 
0 . 000004 
-0 . 000126 
-0 . 000002 
0 . 000000 


-0 . 000005 
-0.003616 
0.000004 
-0,003068 
0.000002 
0.000000 


5.6435  HZ 

0 . 000000 
0 . 000000 
0 . 000000 


7.5861  HZ 

0.000003 

-0.000085 

O.OOOOIB 


.5875  HZ 

-0.000001 

0.000019 

0.000121 


.5877  HZ 

0.000000 

0.000000 

0.000000 


6035  HZ 

0. 000546 
0.004531 
0.005404 


A- 33 


MODE 


20 


FREO 


9.7917  HZ 


MOI'E  = 


MODE  = 


MODE  = 


MODE  = 


0 . 000096 
0.000011 
-0.000001 
0.000015 
-0.000136 
0,000002 


-0.000271 
-0.000036 
-0.000001 
-0.000035 
0 . 000262 
0.000002 


0.003720 

0.010713 

0.011700 


FREO  = 

9.9097  HZ 

-0.347569 
-0.000054 
0.000000 
-0.000077 
0 .497623 
0 . 000000 

0.999941 

0.000154 

0.000000 

0.000157 

-1.018807 

0.000000 

6.283185 

0.000000 

0.000000 

FREQ  = 

9.9102  HZ 

0.000054 
-0,347462 
0 . 000000 
-0.497849 
-0.000077 
0.000000 

-0.000155 

1 .000165 
0.000000 
1.018563 
0.000157 
0.000000 

0.000000 

0.000000 

724.752747 

FREQ  =  10.6229  HZ 

0.000000 
0,000010 
-0.000004 
0.000013 
-0 . 000002 
-  0 . 000000 

-0.000004 
-0.000077 
-0 . 000004 
-0.000061 
0 . 000004 
-  0 . 000000 

0.000001 

0.000000 

0.000000 

-  0 . 000002 
-0.000000 
-0.000001 
-0 , ooooos 

0 . 000014 
0.000004 

0 . 000031 
0 . 000002 
-0 . OOCOOl 
0 . 000005 
-0 . 000058 
0.000004 

0.000000 
0 .000000 
0 . 000003 

A- 34 


Cotiibiried  Model  Risid  Rody 


F  roF-ert  ies 
Mess 

Ceriter  of  nisss 

IneTtia  Matrix 
about  Mass  Center 


9220 

.0  ks 

<  0 

»  0.81345 

»  6.833 

)  ni 

12 

12690 

0 

0 

0  1279220 

-63753 

k  g  ( ni  f  f  2  ) 

0 

-63753 

401329 

A 


-i 
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TABLE  6 


STRENGTH  MODEL  MODES  OF  COMBINED  STRUCTURE  MODEL 
(rero  frepuency  nodes  oniitted). 

All  diselecerients  sre  iri  nieters. 

Modes  ere  riortielised  to  uriity  serierelieed 

The  following  tsble  diseleys  niodel  diselacen.erits.  The 
is  3s  showri. 


ACTUATOR  SENSOR 


LOS  MOTION 


MODE  =  N 


FREQ  =  X.XXXX  HZ 


X  Trerisletiori 

Y  T  r  e  ri  s  1  e  t  i  o  Pi 
Z  TrerislstioPi 
X  Rotation 

Y  Rotatiopi 
Z  Rotation 


MODE  = 


FREQ  =  O.lZOd  HZ 


-0.004151 
--0. 001615 
-0.001004 
0.000061 
-0.002154 
0.000085 


0.002514 

-0.001423 

-0.001004 

0.000211 

-0.002189 

0.000085 


-0.042006 

-0.063311 

0.000625 


MODE 


FREQ 


0.2342  HZ 


■: . ocoooi 

0. CO 00 80 
0.004126 
0 . 000080 
-0 . 000037 
0.000164 


0.000198 

0.000332 

0.004126 

0.000025 

-0.000020 

0.000164 


-0.000001 
-0.000017 
0 . 000733 


MODE  =  3 


FREQ  =  0.3487  HZ 


0.002358 
-0 . 001465 
-0. 002632 
0 . 000682 
0.001155 
-0 . 000989 


0.005213 

0.000650 

-0.002632 

0.001872 

0.002432 

-0.000989 


0 . 004645 
-0 . 347391 
0 . 000003 


ni  3  £  s  . 
f  0  r  ni  s  t 
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MOriE 


A 


FREQ 


0.3816  HZ 


0.000510 
0.000775 
0 . 000456 
-0.000567 
0 . 000346 
-0 . 000005 


0 .001425 
-0.000984 
0 . 000456 
-0.001338 
0.00073? 
-0.000005 


MODE 


FREQ  = 


0 . 000156 
0 . 000404 
0 . 000138 
0.000335 
0.000104 
0.000001 


-0 . 000429 
0.000624 
0.000138 
0.000910 
-0.000226 
-0.000001 


MODE  =  6 


FREQ  = 


-0 . 000031 

-  0  .  C  0  0  2  4  4 

-  0 . 000140 
0 . 000218 
0 . 000328 

-0.000047 


0.004777 
0.000431 
-0.000140 
0.000497 
0.001428 
-0 . 000047 


MODE 


FREQ  = 


-0.001443 
0 . 0034RF 
0.003E  , 

- 0 . 003941 
-0 .001318 
0.001072 


-0 . 004361 
-0.008700 
0 . 003802 
-0.009132 
-0.002693 
0.001072 


MODE  =  8 


FREQ  =  0 


0 .0041  8  5 
0 . 008542 
0 . 005390 
-0 . 01 2399 
0 . 005348 
-0 . 000630 


0.018008 
-0.031320 
0 . 005390 
-0 . 028668 
0 .011262 
-0.000630 


0 . 006990 
0 . 006862 
0.001338 


0.3826  HZ 

-0.002310 

-0.001646 

0.002292 


>.3967  HZ 

-0 . 000429 
0 . 000624 
0.000138 


.4120  HZ 

-0.552454 

0.831537 

0.000625 


.4351  HZ 

0 . 000002 
-0 . 000024 
0.000733 


A- 37 


MOriE  = 

9 

FREQ  = 

0.5001  HZ 

T 

-0 . 06J6ZI 

0 . 0421 60 
-0 . 000373 
-0.1 9304 1 
-0.29171 1 
-0.000021 

-0 . 832899 
-0-551320 
-0 . 000373 
-0.404670 
-0.611541 
-0.000021 

-0.013364 

0.999935 

0.000003 

MOPE  = 

10 

FREQ  = 

0.5002  HZ 

>1 

-0 . 042006 
-0 . 063311 

0 . 000625 
0.291072 
-0. 193555 
-0.000027 

-0 . 552454 
0.831537 
0.000625 
0.610114 
-0.405768 
-0.000027 

0.006990 

0.006862 

0.001338 

k 

MODE  = 

1 1 

FREQ  = 

0.5426  HZ 

0.000219 
-0 . 002005 
0.006337 
-0 . 022774 
-0 . 003232 
-  0.000066 

-0 . 008563 
-0.071448 
0.006337 
-0.047965 
-0.006774 
-0 . 000066 

-0.000045 

-0.000078 

-0.000227 

MODE  = 

1  r 

FREQ  - 

0.6425  HZ 

a 

-0 . 005239 
-0.000241 
0.001141 
0 . 001 626 
0.010139 
-0 . 000562 

0 . 022999 
0.004777 
0.001141 
0 . 002940 
0 . 020957 
-0.000562 

0 . 000000 
0.000000 
0.000000 

h  C'  1 1 E  = 

13 

FREQ  = 

1.4667  HZ 

1 

1 ' 

^9 

0.000020 
-0.014161 
-0 . 001 1 ®2 
0 .019774 
-0 . 00001 4 
O’ .  00000  3 

-0 . 000022 
0 . 055561 
-0.001192 
0.035236 
-0.000021 
0 . 000003 

-0 . 000429 
0 . 000624 
0.000138 

A'  i" 
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fiCOSS  FIFTEEN  (ACTIVE  CONTROL  OF  SPACE  STRUCTURES) (U) 
CONTROL  DVNAHICS  CO  HUNTSVILLE  AL  S  H  SELTZER  ET  AL. 
OCT  82  RADC-TR-82-198  F20682-81-C-0179 
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